AD  A077  689 

^  *  A  ^  ^  ^ 


1.  Report  No. 

NASA-CR  XS9697 


3.  Recipient's  Catalog  No. 


4.  Title  end  Subtitle 


Investigation  of  Critical  Burning  of  Fuel  Droplets  # 


6.  Performing  Organization  Code 


7.  Author(s) 


8.  Performing  Organization  Report  No 


Faeth 


9.  Performing  Organization  Name  end  Addrete 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania  16802 


13.  Type  of  Report  and  Period  Covered 
Final  Report 


National  Aeronautics  and  Space  Administration 
Washington,  D.  C.  20546 


14.  Sponsoring  Agency  Code 


16.  Abstract 

^  A  description  of  the  research  conducted  under  this  grant  is  given  including  theoretical  and 
experimental  methods,  and  major  results.  Three  areas  are  discussed:  (1)  combustion  of 
bipropellant  droplets  with  particular  emphasis  on  high  pressure  conditions  where  the  drop 
approaches  or  exceeds  its  thermodynamic  critical  point  during  combustion;  (2)  combustion  of 
hydrazine  fuels,  treating  hybrid  drop  combustion  in  combustion  gas  environments  and  the 
open-loop  response  of  hydrazine  strands,  drops  and  sprays  to  imposed  pressure  oscillations 
and  (3)  combustion  of  sprays,  where  the  development  of  a  locally  homogeneous  flow  model  of 
spray  combustion  is  described  and  evaluated  with  a  variety  of  single-  and  two -phase  flows,  j 


18.  Distribution  Statement 

Unclassified  -  Unlimited 


Drop  combustion 
High  pressure  combustion 
Spray  combustion 
Hydrazine  fuel  combustion 
Burning  rate  response 


19.  Security  Oanif.  (of  this  report) 


20.  Security  Clatsif.  (of  this  paga) 
Unclassified 


Unclassified 


For  sale  by  the  National  Technical  Information  Service.  Springfield.  Virginia  22161 


NASA-C-168  (Rev.  10-75) 


Final  Report 


Under  NASA  Grant  NGR  39-009-077 


on 

Investigation  of  Critical  Burning  of  Fuel  Droplets 


by 


G.M.  Faeth 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania 


prepared  for 


National  Aeronautics  and  Space  Administration 
Technical  Management 
NASA  Lewis  Research  Center 
Cleveland,  Ohio 

Dr.  Richard  J.  Priem,  Technical  Officer 


Investigation  of  Critical  Burning  of  Fuel  Droplets 


Final  Report 
by 

G.M.  Faeth 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 


SUMMARY 

Three  areas  of  research  were  studied,  as  follows: 

Combustion  of  Bipropellant  Drops.  Observations  were  made  of 
the  combustion  of  drops  in  air,  using  an  apparatus  operating  under 
zero-g  conditions.  High  pressures  were  considered,  where  the  drop 
approaches  or  exceeds  its  thermodynamic  critical  point  during  combus¬ 
tion.  Predictions  of  combustion  lifetime  for  supercritical  combustion 
using  transient  combustion  models,  were  in  good  agreement  with 
measurements.  For  subcritical  combustion,  it  was  found  that  high 
pressure  phenomena  and  the  solubility  of  combustion  product  gases 
in  the  liquid  phase  must  be  considered  if  accurate  predictions  of 
drop  wet-bulb  temperatures  and  the  conditions  required  for  critical 
combustion  are  to  be  obtained.  Heavy  hydrocarbons,  burning  in 
air,  reach  critical  combustion  conditions  at  pressures  on  the  order 
of  twice  the  critical  pressure  of  the  pure  fuel.  Light  hydrocarbons 
tend  to  reach  critical  combustion  conditions  at  pressures  nearer 
their  own  critical  pressure. 

Measurements  of  drop  combustion  rates  in  flowing  air  and  com¬ 
bustion  gases  were  also  completed.  The  pressure  range  of  the  experi¬ 
ments  extended  from  atmospheric  pressure  to  pressures  approaching 
the  critical  combustion  condition.  Predictions  employing  constant 
and  variable  property  combustion  models,  low  and  high  pressure  phase 
equilibrium  models,  and  a  correction  for  convection,  were  compared 
with  these  measurements.  Burning  rates  were  predicted  with  an  average 
error  of  30%  using  average  properties  in  the  convection  correction. 

The  major  difference  between  the  results  of  the  low  and  high  pressure 
phase  equilibrium  models  was  that  the  high  pressure  version  predicts 
a  higher  pressure  for  supercritical  combustion  -  in  better  agreement 
with  the  measurements.  The  use  of  constant  or  variable  property  gas 
phase  models  had  little  impact  on  the  results,  as  long  as  any 
constant  property  was  selected  at  an  arithmetic  average  condition. 

Combustion  of  Hydrazine  Fuels.  The  combustion  of  hydrazine  and 
its  derivatives  as  drops  was  studied  in  combustion  gas  environments 
at  atmospheric  pressure.  When  oxygen  is  present  in  the  gas  around 
the  drop  a  hybrid  flame  structure  is  observed,  with  an  inner  decom¬ 
position  flame,  near  the  surface,  surrounded  by  an  oxidation  zone. 


A  simplified  model  of  this  process  was  developed  which  was  capable 
of  correlating  the  data  of  this  and  other  studies,  with  average 
errors  less  than  20%. 

The  combustion  response  of  hydrazine  to  imposed  pressure  oscil¬ 
lations  was  investigated  both  theoretically  and  experimentally  in 
a  strand  burning  apparatus.  The  theory  was  developed  as  a  perturbation 
analysis  with  the  zero-order  results  yielding  mean  combustion  character 
istics  and  the  first-order  results  yielding  the  linear  combustion 
response.  The  measurements  indicated  linear  response  for  pressure 
oscillation  amplitudes  less  than  15%.  The  analysis  gave  good  pre¬ 
dictions  of  both  steady  combustion  and  combustion  response  character¬ 
istics,  correctly  predicting  both  the  magnitude  and  phase  of  the 
burning  rate  response.  Both  liquid  and  gas  phase  transient  effects 
were  considered  in  the  analysis.  The  results  indicate  that  inter¬ 
actions  between  the  combustion  process  and  liquid  phase  transient 
effects  are  most  important  for  conventional  combustion  processes 
with  hydrazine.  Response  levels  sufficient  to  drive  combustion 
instability  were  found. 

The  strand  results  were  theoretically  extended  to  the  more 
practical  case  of  drops  and  sprays.  It  was  found  that  the  response 
of  monopropellant  hydrazine  drops  and  sprays  can  exceed  unity  for 
drop  sizes,  operating  conditions  and  frequencies  representative  of 
high  frequency  combustion  instability.  The  response  is  greatest  for 
large  drops,  having  mean  temperature  gradients,  at  high  pressures. 
Response  factors  were  less  than  unity  for  drops  at  their  wet -bulb 
state.  The  response  of  a  monodisperse  spray  is  less  than  the  response 
of  the  injected  drop-size  due  to  the  reduced  response  of  small  drops 
near  the  end  of  their  lifetime.  Response  of  polydisperse  sprays 
depends  strongly  on  the  initial  drop  size  distribution;  however,  even 
a  small  percentage  of  large  drops  can  provide  appreciable  response. 

Combustion  of  Sprays.  A  spray  combustion  model  was  developed  and 
systematically  evaluated  employing  a  sequence  of  experiments  as  follows 
constant  density  jet,  variable  density  jet,  evaporating  spray,  com¬ 
busting  jet  and  combusting  spray.  The  experimental  conditions  involved 
jets  or  two-phase  flows  in  a  stagnant  environment.  Measurements  at 
atmospheric  pressure  consisted  of  mean  and  turbulent  velocities, 
using  the  laser  Doppler  anemometer;  mean  temperatures,  using  fine 
wire  thermocouples;  mean  composition,  using  isokinetic  sampling  and 
analysis  by  gas  chromatography;  and  drop  size  distributions,  using  a 
particle  impaction  method.  Data  at  elevated  pressures  consisted  of 
photographs  of  spray  and  flame  boundaries.  Test  pressures  extended  to 
supercritical  conditions. 

The  theoretical  model  was  based  on  the  assumption  of  locally  homo¬ 
geneous  flow  (LHF)  for  the  two-phase  flows.  This  assumption  implies 
that  local  gas  and  liquid  velocities  are  the  same  and  that  the  phases 
are  in  thermodynamic  equilibrium,  at  each  point  in  the  flow.  Flow 
characteristics  were  based  on  a  k-e-g  turbulence  model.  Under  these 
assumptions,  the  various  flows  only  differ  by  their  equation  of 
state,  which  relates  local  mixture  fraction  to  other  properties  in  the 
flow,  e.g.  density,  temperature,  composition,  etc. 


iii 


Using  a  well-established  set  of  turbulence  model  constants,  the 
model  predictions  were  in  excellent  agreement  with  the  measurements 
for  all  single-phase  flows,  both  noncombusting  and  combusting.  The 
model  provided  a  useful  qualitative  picture  of  the  flow  in  evaporating 
and  combusting  sprays  at  atmospheric  pressure  (maximum  errors  in  mean 
velocity  and  mean  mixture  fraction  were  in  the  range  20%)  for  sprays 
having  a  Sauter  mean  diameter  on  the  order  of  30  urn.  The  errors 
result  from  effects  of  slip  and  loss  of  thermodynamic  equilibrium 
between  the  phases.  Analysis  of  drop-life  histories  suggested  that  a 
Sauter  mean  diameter  less  than  10  ym  would  be  required  for  quantitative 
accuracy  at  the  atmospheric  pressure  test  conditions.  The  comparison 
between  theory  and  experiment  was  improved  at  high  pressures,  although 
additional  experimental  evaluation  is  needed  in  order  to  obtain  a 
more  definitive  evaluation. 

The  advantage  of  the  LHF  model  is  that  detailed  information  on 
spray  conditions  are  not  needed  for  a  prediction.  Effects  of  slip 
and  loss  of  thermodynamic  equilibrium  cause  the  LHF  model  to  overestimate 
the  rate  of  flow  development,  thus  the  model  represents  a  lower  bound 
for  all  possible  sprays.  Given  the  predictions  of  the  LHF  model,  the 
potential  importance  of  slip,  etc.  can  be  assessed  by  completing  drop- 
life  history  calculations,  using  the  LHF  results  to  specify  the  local 
environment  of  drops  in  the  spray. 
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NOMENCLATURE 

Description 

gravitational  acceleration 
Damkohler  number 
constants  in  turbulence  model 
injector  passage  diameter 
drop  diameter 

activation  energy,  dimensionless  activation  energy 
mixture  fraction 

square  of  mixture  fraction  fluctuations 

turbulence  kinetic  energy 

burning  rate  constant 

injector  mass  flow  rate 

jet  momentum 

reaction  order 

pressure 

reduced  pressure,  drop  response  factor 
spray  response  factor 

mass  of  oxidant  per  unit  mass  of  fuel  for  stoichiometric 
combustion,  radial  distance 

time 

temperature 
axial  velocity 
radial  velocity 

effective  radial  velocity,  Eq.  (6) 

axial  distance 

mass  fraction  of  species  i 

dissipation  rate  of  turbulent  kinetic  energy 


yt  turbulent  viscosity 

v  stoichiometry  parameter 

p  density 

o  turbulent  Prandtl/ Schmidt  number 

X  transport  parameter 

co  dimensionless  frequency 


Subscripts 

A 
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Ff 

Fg 
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o 
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air 

centerline 
injected  fluid 

injected  fluid  in  liquid  phase 
injected  fluid  in  gas  phase 
initial  state 
initial  state 
oxygen 

ambient  state 
steady  flame  condition 


Superscripts 


mean  quantity 
fluctuating  quantity 
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I.  INTRODUCTION 

This  investigation  considers  aspects  of  the  general  problem  of 
spray  combustion .  The  research  was  primarily  motivated  by  a  need  for 
better  understanding  of  stable  and  unstable  combustion  in  liquid-fueled 
rocket  engines.  However,  the  results  also  have  application  to  other 
liquid-fueled  combustion  systems,  e.g.,  gas  turbines,  furnaces  and 
Diesel  engines. 

Three  areas  of  research  were  studied  under  the  grant: 

1.  Combustion  of  bipropellant  drops.  In  this  area,  particular 
emphasis  was  placed  on  nigh  pressure  combustion  of  individual  drops, 
under  conditions  where  the  liquid  drop  approaches  or  exceeds  its 
thermodynamic  critical  point.  Theory  was  developed  to  assist  the 
interpretation  of  experimental  results  and  to  establish  methods  for 
computing  drop  behavior  for  conditions  that  were  not  explicitly 
tested.  The  measurements  involve  pure  hydrocarbon  and  alcohol  fuels,  in 
air  and  combustion  gas  environments,  at  pressures  up  to  13.8  MPa 

(2000  psia) . 

2.  Combustion  of  hydrazine  fuels.  Combustion  of  hydrazine,  MMH  and 
UDMH  as  drops  was  studied  both  theoretically  and  experimentally.  This 
involved  basic  measurements  of  combustion  rates  of  drops  in  combustion 
gas  environments,  and  modeling  the  process  to  include  both  decomposition 
and  oxidation  phenomena.  The  combustion  response  of  liquid  hydrazine 

to  imposed  pressure  oscillations  was  also  studied.  Predicted  and 
measured  response  for  burning  hydrazine  strands  were  .compared.  Based 
on  these  results,  predictions  of  drop  and  spray  response  were  obtained. 

3.  Combustion  of  sprays.  The  final  stages  of  the  investigation  con- 
sidered  the  development  of  a  model  for  spray  combustion  and  obtaining 
experimental  results  to  evaluate  this  model.  High  pressure  conditions 
where  phenomena  associated  with  the  liquid  approaching  its  thermo¬ 
dynamic  critical  point  are  significant  was  considered  as  well  as  low 
pressure  operation.  The  model  was  evaluated  systematically  considering: 
constant  density  jets,  variable  density  jets,  combusting  jets,  evapor¬ 
ating  sprays,  and  combusting  sprays. 

In  this  report,  each  area  of  the  research  is  briefly  described, 
and  the  major  findings  are  discussed.  Further  details  may  be  found  in 
specific  reports  (1-18)*,  papers  (19-38)  and  theses  (39-49)  resulting 
from  this  grant. 

II.  COMBUSTION  OF  BIPROPELLANT  DROPS 
2.1  Introduction 


A  central  problem  that  mu'- 1  be  resolved  prior  to  any  rational 
treatment  of  spray  combustion  involves  the  transport  characteristics 
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Numbers  in  parenthesis  denote  references. 
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of  individual  drops  within  a  spray.  This  phase  of  the  investigation 
was  concerned  with  processes  of  heat  and  mass  transfer  to  individual 
drops.  A  variety  of  test  conditions  were  examined;  however,  the  main  ’ 

emphasis  was  placed  on  hydrocarbon  and  alcohol  fuels,  in  combustion 
gas  environments  at  high  pressures  (where  operation  near  the  thermo¬ 
dynamic  critical  point  of  the  fuel  must  be  considered) . 

Two  topics  were  considered  in  this  portion  of  the  research, 
as  follows: 

1 •  Supercritical  drop  combustion.  This  topic  involves  conditions 
where  the  drop  exceeds  its  thermodynamic  critical  point,  and  burns  in 
a  transient  manner,  similar  to  a  puff  of  gas.  Experiments  were  con¬ 
ducted  to  observe  drop  combustion  under  these  conditions.  The  measure¬ 
ments  were  compared  with  theoretical  results  in  the  literature. 

2.  High  pressure  drop  combustion.  This  area  involved  the  determination 
of  the  pressure  levels  required  for  supercritical  combustion,  and  heat 
and  mass  transfer  rates  in  the  near  critical  pressure  regime.  Theory 
was  developed  for  these  parameters,  and  compared  with  the  measurements. 

The  main  results  obtained  in  each  of  these  areas  is  described  in 
the  following  sections.  « 

2.2  Supercritical  Drop  Combustion 

Introduction.  Theoretical  investigations  by  Spalding  (50), 

Wieber  (51),  Brzustowski  (52)  and  Rosner  (53)  indicated  that  at 
sufficiently  high  pressures,  a  fuel  drop  can  exceed  its  thermodynamic 
critical  point  during  combustion.  When  this  occurs,  the  quasisteady 
drop  combustion  theories ,  employed  with  good  success  to  predict 
liquid- fueled  rocket  engine  performance  at  moderate  pressures  (54),  no 
longer  represent  the  combustion  process  and  a  new  burning  regime  must 
appear.  The  pressure  regime  for  this  behavior  is  in  the  range  of 
operation  of  some  liquid-fueled  rocket  engines  and  Diesel  engines. 

Some  modem,  high  pressure  ratio  gas  turbine  combustors  approach  this 
region  of  operation  as  well. 

Upon  reaching  its  thermodynamic  critical  point,  a  droplet  becomes 
a  puff  of  gas.  Spalding  (50)  theoretically  considered  combustion  in 
this  case  by  approximating  the  droplet  vapor  as  an  instantaneous  point 
source  of  fuel.  Laminar  combustion  and  spherical  symmetry  were  assumed 
in  the  analysis.  The  flame  was  taken  to  be  infinitely  thin,  with  fuel 
and  oxident  combining  in  stoichiometric  proportions.  Physical 
properties  were  assumed  to  be  constant.  Rosner  (53),  modified 
Spalding's  analysis  to  account  for  the  initially  finite  dimensions  of 
the  fuel  vapor  pocket. 

P 

Both  these  analyses,  and  the  theories  used  to  suggest  the  onset 
of  critical  combustion,  contain  many  approximations.  Therefore,  the 
main  objective  of  the  present  investigation  was  to  obtain  experimental 
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results  to  test  the  predictions.  Only  a  limited  description  of  the 
investigation  is  presented  here,  further  details  can  be  found  in 
Refs.  1-3,  19-21,  39  and  41. 

Apparatus .  A  sketch  of  the  test  apparatus  appears  in  Fig.  1 
The  apparatus  consisted  of  a  pressure  vessel  fitted  with  windows  to 
allow  motion  pictures  of  the  combustion  process.  A  drop  was  supported 
at  the  center  of  the  chamber  on  either  a  quartz  fiber  or  a  thermo¬ 
couple  junction.  The  drop  was  ignited  with  a  hot  wire.  The  gaseous 
environment  of  the  drop  consisted  of  various  mixtures  of  oxygen  and 
nitrogen . 

As  the  critical  point  of  a  drop  is  approached,  surface  tension 
approaches  zero,  and  the  supported  drop  technique  is  not  feasible 
since  the  drop  would  fall  from  the  probe.  This  problem  was  avoided 
by  conducting  the  experiment  in  a  freely  falling  apparatus,  which 
provided  a  test  time  of  approximately  Is  under  nearly  zero-gravity 
conditions.  Data  was  obtained  on  drop  temperature  and  flame  shape 
variations  during  combustion,  as  well  the  time  for  complete 
combustion. 

Results  and  Discussion.  Figure  2  is  an  illustration  of  drop 
temperature  variation  as  a  function  of  time  after  ignition  at  various 
pressures.  The  inflection  of  the  temperature  trace  a  short  time  after 
ignition  corresponds  to  the  point  when  the  thermal  disturbance  due 
to  the  presence  of  a  flame  propagates  to  the  thermocouple  position. 

At  lower  pressures,  curve  A  of  Fig.  2,  the  temperature  of  the 
drop  levels  out  at  the  wet  bulb  temperature  (the  condition  where  all 
the  heat  transferred  from  the  flame  to  the  drop  is  utilized  for  the 
heat  of  vaporization  of  the  evaporating  fuel).  The  temperature  then 
rises  again  after  all  the  liquid  is  gone.  This  behavior  is  charac¬ 
teristic  of  subcritical  combustion,  treated  by  the  vaporization 
controlled  models (54) . 

As  the  pressure  is  increased,  curve  B  of  Fig.  2,  the  time  spent 
at  the  wet  bulb  state  decreases.  Eventually,  a  condition  is  reached 
where  the  drop  temperature  increases  continuously,  curve  C.  This  is 
indicative  of  supercritical  combustion,  which  persists  for  all  higher 
pressures.  In  this  case,  the  fuel  becomes  gradually  more  vapor- like 
as  its  temperature  increases  and  the  bulk  of  the  liquid  does  not 
evaporate  from  a  liquid  surface  as  it  does  at  low  pressures. 

Tests  were  conducted  undeT  supercritical  combustion  conditions  and 
measurements  were  made  of  the  time  of  combustion,  for  comparison  with 
the  theories  of  Refs.  50  and  S3.  A  portion  of  these  results  is  illus¬ 
trated  in  Fig.  3.  Combustion  lifetime,  corrected  to  a  fixed  initial 
drop  size,  is  plotted  as  a  function  of  a  stoichiometric  factor. 
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Fig,  1  Sketch  of  the  zero-gravity  test  chamber. 
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Droplet  temperature  measurements  for  n-decane  in  air  at  various 
pressures.  Pressures  and  initial  drop  diameters:  Curve  A,  0.79 
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where  is  the  ambient  oxident  concentration  and  r  is  the  mass  of 
oxident  required  to  consume  a  unit  mass  of  fuel.  The  difference  be¬ 
tween  the  two  theories  is  not  large  for  the  test  range  of  Fig.  3.  The 
experimental  results  are  in  good  agreement  with  predictions. 

Other  factors,  such  as  initial  drop  size,  flame  shape  and  total  < 

pressure,  were  also  examined  during  the  study.  The  predictions  were 
satisfactory  with  respect  to  these  variables  as  well. 

Conclusions .  This  study  provided  the  first  experimental  obser¬ 
vations  of  supercritical  combustion  of  individual  drops.  No  unusual 
change  in  the  structure  of  the  combustion  zone  was  observed  as  the 
drop  passed  through  its  critical  condition.  In  general,  predictions  of 
the  two  theoretical  models  were  in  reasonably  good  agreement  with  the 
observations . 

2.3  High  Pressure  Drop  Combustion 
2.3.1  Zero-Gravity  Combustion 

Introduction .  With  supercritical  combustion  resolved, 
at  least  in  the  absence  of  forced  convection,  attention  turned  to  de¬ 
termining  the  conditions  required  for  supercritical  combustion  to  t 

appear.  Wieber  (51)  had  earlier  computed  total  pressures  required  for 
supercritical  evaporation  of  oxygen  and  n-heptane  drops  in  a  high 
temperature  inert  gas.  However,  this  analysis  did  not  allow  for  the 
influence  of  high  pressures  on  the  thermodynamic  properties  of  the 
fuels,  and  particularly,  the  increased  solubility  of  liquid  fuels  for 
combustion  gases  at  high  pressures.  Manrique  and  Borman  (55)  suggested 
that  these  effects  have  a  significant  influence  on  drop  equilibrium 
conditions  at  high  pressures. 

In  the  light  of  these  findings,  high  pressure  drop  combustion  was 
considered  both  theoretically  and  experimentally.  Emphasis  was  placed 
on  determining  wet-bulb  temperatures  during  subcritical  drop  combustion. 

Complete  details  of  this  portion  of  the  study  can  be  found  in  Refs. 

4-6,  21-24,  and  41. 

Apparatus .  The  apparatus  was  essentially  the  same  as  described  in 
Section  2.2.  Drop  temperatures  were  measured  with  thermocouples  con¬ 
structed  of  25,  50  and  75  pm  wires.  Results  with  the  smallest  wires 
suggested  that  errors  were  small  using  50  pm  wires  and  these  were 
employed  for  the  bulk  of  the  tests. 

Theory.  A  theoretical  model  for  the  combustion  process  at  high 
pressures  was  developed.  There  are  two  aspects  to  the  model,  a  gas 
phase  combustion  model  and  a  model  for  phase  equilibrium  at  the  drop 
surface.  9 

The  gas  phase  combustion  model  was  an  extension  of  the  quasisteady 
drop  burning  model  of  Goldsmith  and  Penner  (56) .  This  model  assumes 
spherical  symmetry,  quasisteady  combustion  at  the  wet-bulb  state  and 
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an  infinitely  thin  diffusion  flame.  The  model  provides  for  variable 
properties  in  the  gas  phase.  The  extension  completed  in  this  in¬ 
vestigation  involved  consideration  of  the  effect  of  dissolved  gas 
evaporation,  separate  determination  of  the  concentrations  of  various 
gas  phases  species,  and  allowing  for  variable  specific  heats  of  all 
species . 

Two  models  were  employed  for  computing  phase  equilibrium 
conditions  at  the  surface  of  the  drop.  The  simplest  model  corres¬ 
ponded  to  the  usual  treatment  of  phase  equilibrium,  valid  only  at 
low  pressures  (54) .  The  second  model  provided  for  solubility  and  high 
pressure  effects,  using  a  modified  Redlich-Kwong  equation  of  state 
developed  by  Prausnitz  and  Chueh  (57) . 

Results  and  Discussion.  Figure  4  is  an  illustration  of  predicted 
and  measured  wet -bulb  temperatures  as  a  function  of  pressure  for  the 
combustion  of  n-octane  drops  in  air.  The  boiling  curve,  and  the  pre¬ 
dictions  from  the  low  and  high  pressure  phase  equilibrium  model  are 
shown  on  the  figure.  Uncertainty  bands  are  indicated  on  the  theoretical 
results.  Greatest  uncertainties  were  due  to  a  transport  parameter  xF* 
which  is  not  known  very  accurately.  The  theoretical  curves  are 
terminated  where  steady  burning  solutions  no  longer  exist.  The  experi¬ 
mental  results  are  terminated  when  an  inflection  could  no  longer  be 
observed  in  the  plot  of  drop  temperature  variation  with  time,  c.f. 

Fig.  2.  The  supercritical  region  occurs  at  pressures  higher  than  this, 
and  follows  the  unsteady  combustion  mechanism  described  in  Section  2.2. 

The  high  pressure  model  is  in  reasonably  good  agreement  with  the 
measurements  in  Fig.  4,  both  with  respect  to  the  wet-bulb  temperature 
and  the  pressure  range  required  for  supercritical  combustion.  The  low 
pressure  theory  underestimates  the  pressure  required  for  critical 
combustion  by  roughly  a  factor  of  2. 

Measured  and  predicted  pressures  for  supercritical  combustion  are 
summarized  in  Table  1  for  n-pentane,  n-octane  and  n-decane.  The  high 
pressure  theory  provides  the  best  prediction  of  supercritical  com¬ 
bustion  pressure.  It  is  also  seen  that  the  critical  combustion  pressure 
tends  to  approach  the  critical  pressure  of  the  pure  fuel  as  the 
molecular  weight  of  the  fuel  decreases. 

The  general  characteristics  of  steady  evaporation  at  high  pressures 
are  illustrated  in  Fig.  5  The  plot  illustrates  the  locus  of  steady 
evaporation  states  from  the  gas  phase  model  and  the  locus  of  phase 
equilibrium  states  from  the  phase  equilibrium  theory.  The  two  loci 
yield  fuel  mole  fractions  in  the  gas  phase  at  the  liquid  surface,  as 
a  function  of  liquid  temperature,  with  reduced  pressure  as  a  parameter. 
As  a  droplet  heats-up,  after  ignition,  it  follows  the  locus  of  phase 
equilibrium  states.  At  the  point  where  the  two  loci  intersect,  drop 
heating  ceases,  and  the  steady  evaporation  period  of  the  drop  begins. 

As  the  pressure  is  increased,  a  condition  is  reached  where  the  two  loci 
are  just  tangent.  This  is  the  highest  pressure  where  steady  evaporation 
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Theoretical  and  experimental  steady-burning  temperatures  for  n-octane 
droplets  burning  in  air  under  zero-gravity  conditions. 
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Fig.  5  Predictions  of  droplet-surface  fuel  mole  fractions  as  a  function 

of  liquid  temperature  for  n-decane  evaporating  in  nitrogen  at  3000  K. 
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Table  1 

Measured  and  Predicted  Pressures  Required  for  Supercritical  Combustion 


Fuel 

n-Pentane 

n-Octane 

n-Decane 

Measured  critical  burning 
pressure  (atm) 

- 

58-61 

44-47 

Measured  critical  P 

r 

- 

2. 4-2.5 

2. 1-2. 3 

Predicted  critical  Pr,  low 
pressure  theory*3 

1.1-1. 2 

1.2-1. 4 

1.2-1. 5 

Predicted  critical  P  ,  high 
pressure  theory*3  r 

1.2-1. 7 

1.7-2. 3 

1.8-2. 8 

Pr  *  (critical  burning  pressure)/ (thermodynamic  critical 
pressure  of  pure  fuel) 

variation  due  to  +_  20%  variation  in  xF 


r 
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is  observed.  At  this  limit,  the  drop  evaporates  at  a  thermodynamically 
subcritical  state.  At  higher  pressures,  heat  input  to  the  drop  from 
the  flame  always  exceeds  the  energy  required  to  vaporize  the  fuel,  and 
drop  heating  continues  along  the  locus  of  phase  equilibrium  states. 

When  the  locus  of  critical  mixing  states  is  crossed,  the  material 
become  gas-like  and  subsequently  burning  is  unsteady  supercritical 
combustion. 

Conclusions .  The  results  indicate  that  high  pressure  effects  and 
the  solubility  of  combustion  gases  in  the  liquid  phase  must  be  con¬ 
sidered  if  accurate  predictions  of  drop  wet-bulb  temperatures  and  the 
pressure  required  for  critical  combustion  are  to  be  obtained.  For 
heavier  hydrocarbons,  burning  in  air,  the  critical  combustion  pressure 
is  on  the  order  of  twice  the  critical  pressure  of  the  pure  fuel. 

Lighter  hydrocarbons,  however,  tend  to  reach  critical  combustion 
conditions  at  pressures  nearer  their  thermodynamic  critical  pressure. 
Conditions  were  also  found,  at  high  pressures,  where  water  vapor  pro¬ 
duced  by  combustion  condenses  on  the  liquid  surface  during  the  early 
stages  of  drop  heat-up. 

2.3.2  Combustion  under  Convection  Conditions 

Introduction.  The  zero-gravity  experiments  were  useful 
for  examining  critical  combustion  conditions;  however,  attempts  to  ob¬ 
tain  information  on  drop  gasification  rates  were  not  particularly 
successful.  The  heated  zone  around  the  droplet  caused  excessive  bending 
of  background  light  beams;  therefore,  conventional  shadowgraphs  could 
not  be  employed  to  find  the  variation  of  drop  size  with  time,  at  high 
pressures.  Under  zero-gravity  conditions,  the  drop  combustion  process, 
even  at  subcritical  combustion  conditions,  is  intrinsically  unsteady 
(34).  Unsteadiness  increases  the  difficulties  of  data  interpretation. 
Finally,  drop  motion  is  a  significant  effect  for  drops  in  sprays  and 
information  in  this  area  is  needed. 

With  these  problems  in  mind,  a  series  of  experiments  were  conducted 
in  order  to  examine  both  high  pressure  phenomena  and  convection.  As  a 
baseline,  tests  were  conducted  at  atmospheric  pressure,  with  drops 
supported  in  combustion  products  of  a  flat  flame  burner.  Subsequent  tests 
employed  the  porous  sphere  technique,  which  allows  unsteadiness  to  be 
eliminated.  These  results  provided  burning  rates  and  liquid  surface 
temperatures  at  pressures  extending  to  the  supercritical  combustion 
condition.  Both  air  and  combustion  gases  were  employed  as  the  environ¬ 
ment  foT  porous  sphere  combustion.  The  experiments  were  compared  with 
predictions,  using  models  similar  to  those  discussed  in  Section  2.3.1. 

Complete  details  of  this  portion  of  the  investigation  can  be  found 
in  Refs.  8-11,  25,  28,  30,  41,  and  44. 

Apparatus .  Two  different  test  arrangements  were  employed  for  the 
convection  experiments.  The  supported  droplet  technique  was  employed 
for  baseline  tests  at  atmospheric  pressure.  The  test  arrangement  in 
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this  case  is  illustrated  in  Fig.  6.  The  drops  were  supported  on  a  quartz 
filament.  Combustion  gas  environment  was  provided  by  a  flat  flame  burner. 
Various  burner  mixture  ratios  were  employed  to  achieve  a  range  of  tem¬ 
peratures  and  oxident  concentrations  for  the  drop  environment.  The 
burner  could  be  rapidly  moved  under  the  droplet,  using  a  pneumatic 
cylinder,  in  order  to  initiate  the  combustion  process.  The  variation 
of  drop  diameter  with  time  was  obtained  from  motion  picture  shadow¬ 
graphs  . 

The  second  test  arrangement  employed  the  porous  sphere  technique. 

A  sketch  of  this  apparatus  appears  in  Fig.  7.  The  tests  were  conducted 
within  a  water-cooled  pressure  vessel,  fitted  with  windows  to  allow 
observation  of  the  combustion  process.  Droplets  were  simulated  by 
porous  spheres  of  various  sizes.  The  fuel  flow  rate  to  the  sphere 
could  be  varied,  using  a  variable  displacement  pump.  The  steady  com¬ 
bustion  rate  of  the  sphere  was  determined  as  the  point  where  the 
sphere  neither  dripped  liquid  or  dried-off.  A  system  of  burettes  at 
the  pump  inlet  was  used  to  determine  the  steady  burning  rate.  Com¬ 
bustion  gas  environments  were  provided  by  a  burner  housed  at  the  bottom 
of  the  pressure  vessel.  Control  of  the  ambient  conditions  of  the 
sphere  was  similar  to  method  employed  for  the  flat  flame  burner  appa¬ 
ratus.  Tests  were  also  conducted  with  the  sphere  in  an  air  environment 
•  with  this  apparatus.  In  this  case,  a  slow  flow  of  air  was  maintained 

through  the  chamber  in  order  to  prevent  vitiation  of  the  air.  Drop 
surface  temperatures  were  measured,  using  76  um  thermocouples  imbedded 
flush  with  the  surface  of  the  porous  sphere. 

Theory.  The  combustion  model  discussed  in  Section  2.3.1  was  used 
for  comparison  with  the  present  measurements.  Spherical  symmetry  and 
no  convection  was  assumed  for  the  original  model.  Therefore,  the  effect 
of  convection  was  included  as  correction  factor,  multiplying  the  pre¬ 
diction  for  no  convection.  This  is  a  standard  procedure,  widely  used 
in  the  mass  transfer  literature  (34,  51,  54).  The  correlation  used  in 
the  present  work  has  good  asymptotic  properties,  and  properly  agrees 
with  theory  at  the  low  Reynolds  number  limit  (34) . 

As  in  earlier  work  (Section  2.3.1)  both  the  low  and  high  pressure 
phase  equilibrium  models  were  employed  for  the  predictions.  For 
porous  sphere  combustion,  however,  the  boundary  condition  at  the  liquid 
surface  must  be  modified  from  that  employed  for  drop  combustion.  This 
is  due  to  the  fact  that  the  liquid  entering  the  sphere  must  always  be 
heated  to  the  surface  temperature,  while  bulk  heating  of  a  drop  ceases 
during  the  steady  evaporation  period. 

Results  and  Discussion.  Some  typical  results  of  the  baseline  sup¬ 
ported  drop  tests  are  illustrated  in  Fig.  8.  Over  a  limited  diameter 
range,  with  the  drop  at  its  wet-bulb  state,  the  square  of  the  drop 
diameter  varies  linearly  (34) 
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ARC  BACKGROUND  LIGHT  QUARTZ  FILAMENT 


Fig.  6  Sketch  of  the  flat  flame  burner  apparatus. 
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where  K  is  called  the  burning  rate  constant.  The  burning  rate  constant 
is  illustrated  in  Fig.  8  as  a  function  of  fuel  type  for  various  normal 
paraffins.  Ambient  oxygen  concentration  is  a  parameter  on  this  plot. 

The  temperature  level  of  the  tests  is  2530  K.  The  convection  correction 
for  these  conditions  is  less  than  23  percent  of  the  no-convection 
prediction  (25) . 

There  is  fair  agreement  between  theory  and  experiment  in  Fig.  8. 

The  effect  of  oxident  concentration  is  well  represented.  The  effect  of 
fuel  type  is  somewhat  less  satisfactory.  Other  investigators  have  also 
observed  that  conventional  drop  burning  theories  tend  to  overestimate 
measurements  for  heavier  hydrocarbons  (58) .  This  behavior  is  still 
not  fully  resolved,  effects  of  variable  properties  and  gas  phase  fuel 
decomposition  have  been  suggested  as  possible  reasons  for  the  diffi¬ 
culty.  Results  for  alcohols  were  similar  to  the  findings  for  paraffin 
hydrocarbons.  In  general,  the  theory  provides  a  reasonably  sound  means 
of  estimating  the  effects  of  ambient  composition,  temperature  and 
pressure  on  drop  burning  rates  for  a  given  fuel.  Fuel  type  effects, 
however,  are  best  handled  by  matching  the  prediction  with  measurements 
for  at  least  one  condition..  There  is  sufficient  uncertainty  in 
physical  properties  for  the  gas  phase,  even  for  so-called  variable 
property  models,  that  a  proper  match  can  generally  be  achieved  by 
*  adjusting  properties  within  reasonable  bounds. 

Figure  9  is  an  illustration  of  measured  fuel  burning  rates  for 
,  porous  spheres,  burning  in  air,  as  a  function  of  pressure.  Various 

noxmal  paraffins  are  considered,  under  natural  convection  conditions. 
Theoretical  results  are  compared  with  the  measurements,  considering 
both  the  low  and  high  pressure  phase  equilibrium  theories.  The  theo¬ 
retical  curves  are  terminated  at  their  high  pressure  end  when  the 
critical  burning  pressure  is  reached.  The  experimental  results  are 
terminated  at  high  pressures  due  to  problems  with  carbon  spots  foraing 
on  the  surface  of  the  sphere.  Both  theories  are  in  reasonably  good 
agreement  with  the  measurements.  The  major  difference  between  the 
theories  is  that  the  low  pressure  theory  predicts  a  lower  critical 
combustion  condition,  as  noted  earlier  for  drops. 

Results  for  alcohol  fuels  were  similar  to  the  findings  illustrated 
in  Fig.  9.  Since  lower  alcohols  produce  less  soot,  it  was  possible 
to  observe  supercritical  combustion  with  methanol  and  ethanol  at 
pressures  in  the  range  80-100  atm.  For  supercritical  combustion 
conditions,  the  fuel  flow  rate  could  be  varied  arbitrarily  and  no 
dripping  of  fuel  was  observed.  Increasing  the  fuel  flow  rate  simply 
caused  the  flame  to  move  away  from  the  surface.  Similar  effects  are 
observed  for  the  combustion  of  gases  from  porous  spheres. 

The  effect  of  sphere  diameter  was  examined  for  diameters  in  the 
range  6.4-19  mm.  The  theories  accurately  predicted  the  influence  of 
diameter  on  the  burning  rate,  under  natural  convection  conditions. 

The  effect  of  combustion  gas  environments  was  also  considered.  The 
results  were  similar  to  Fig.  9.  In  general,  the  discrepancy  between 
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Fig.  9  Porous  sphere  burning  rates  of  paraffins  in  air,  for  fuel 
inlet  and  ambient  temperatures  of  300  K.  Sphere  diameter. 
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theory  and  predictions  was  greatest  when  the  ambient  temperature  was 
high  and  the  ambient  oxygen  concentration  was  low. 

Figure  10  is  a  comparison  of  the  overall  agreement  between  theory 
and  experiment  for  combustion  gas  environments.  The  data  includes  both 
the  baseline  results  at  atmospheric  pressure  (25) ,  and  high  pressure 
results  (30) .  The  convection  correction  was  substantial  for  high 
pressure  conditions.  The  data  points  shown  on  the  figure  correspond 
to  the  use  of  ambient  properties  in  the  convection  correction,  which 
was  a  procedure  advocated  by  some  investigators  (25) .  Better  agree¬ 
ment  was  obtained  using  average  properties  in  this  correction.  The 
individual  points  when  average  properties  are  used  are  not  shown  on 
Fig.  10,  in  order  to  avoid  cluttering  the  figure;  however,  the  bounds 
of  the  data  are  indicated  by  the  lines  A-A.  The  choice  of  average  or 
ambient  conditions  in  the  convection  correction  for  the  Faeth  and  Lazar 
(25)  data  is  not  very  significant,  since  the  correction  was  small  for 
these  test  conditions.  The  test  conditions  for  all  the  data  illus¬ 
trated  in  Fig.  10  are  summarized  in  Table  2. 

Figure  11  is  an  illustration  of  predicted  liquid  and  gas  phase 
compositions,  at  the  liquid  surface,  for  porous  sphere  combustion  in 
a  combustion  gas  environment.  The  gas  phase  composition  remains  rela¬ 
tively  constant  as  the  total  pressure  is  increased.  However,  the  liquid 
phase  composition  of  dissolved  gases  increases  steadily  with  increasing 
pressure.  The  thermodynamic  critical  point  is  indicated  by  the  condition 
where  the  composition  of  both  phases  is  the  same.  The  results  are  similar 
for  drop  combustion,  except  that  the  highest  pressure  for  subcritical 
combustion  is  not  the  mixture  critical  point,  and  the  composition  of 
both  phases  are  not  the  same  at  this  condition. 

Conclusions.  The  low  and  high  pressure  theories  yield  roughly  the 
same  burning  rate  prediction.  The  majoT  difference  between  these  models 
is  that  the  high  pressure  theory  predicts  a  higher  pressure  for  critical 
combustion.  Predictions  of  critical  combustion  using  the  high  pressure 
model  are  in  reasonably  good  agreement  with  measurements.  The  use  of 
constant  or  variable  properties  in  the  gas  phase  models  had  little 
impact  on  the  results,  as  long  as  any  constant  property  was  selected 
at  an  arithmetic  average  condition.  This  is  particularly  important  for 
convection  corrections,  where  average  properties  should  be  used,  c.f. 

Fig.  10.  Following  these  prescriptions,  burning  rates  for  combustion 
gas  environments  could  be  predicted  with  an  average  error  of  30%,  over 
the  range  of  variables  summarized  in  Table  2.  More  accurate  results 
can  be  obtained  by  matching  theory  and  experiment  for  at  least  one 
operating  condition,  since  existing  theories  are  least  effective  con¬ 
cerning  effects  of  fuel  type. 

Similar  to  the  results  discussed  in  Section  2.3.1,  the  high  pressure 
model  provided  the  best  prediction  of  liquid  surface  temperatures.  Dis¬ 
solved  gas  concentrations  become  very  significant  at  high  pressures, 
e.g. ,  reaching  values  approaching  50  percent  in  Fig.  11.  This  tends  to 
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SYMBOL 


FUEL 
METHANOL 
ETHANOL 
I- PROPANOL 
I- BUTANOL 
l-DECANOL 
n  -PENTANE 
n  -HEPTANE 
i-  OCTANE 
n  -DECANE 
n -TRIDECANE 
n -HEXADECANE 
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Overall  comparison  of  theoretical  and  experimental  gasification 
rates  in  combustion  gas  environments.  Lines  A-A  are  data  bounds 
when  average  properties  are  used  in  the  convection  correction. 
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Predicted  liquid-surface  compositions  for  porous  sphere 
combustion.  Ambient  temperature  1145  K,  Fuel  inlet 
temperature  300  K. 
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Table  2 

Summary  of  Test  Conditions  for  Drop  Combustion 
Rates  in  Combustion  Gas  Environments 


Canada  and  Faeth  (30) 

Faeth  and  Lazar  (25) 

Test  Method 

Porous  Sphere 

Suspended  Drop 

Ambient  Temperature  (K) 

600-1500 

1660-2530 

Ambient  Oxygen  Cone. 

(mole  %) 

0-13.5 

0-41.8 

Pressure  (atm) 

1-40 

1 

Diameter  (pm) 

9500 

1100 

Reynolds  number 

(approach  properties) 

40-680 

1.3-2. 2 
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raise  the  critical  pressure  of  the  mixture  at  the  liquid  surface,  and 
reduce  the  critical  temperature,  from  that  of  most  hydrocarbon  fuels. 
The  net  result  is  that  critical  combustion  occurs  at  pressures  in 
excess  of  the  critical  pressure  of  the  pure  fuel. 

III.  COMBUSTION  OF  HYDRAZINE  FUELS 

3 . 1  Introduction 


Hydrazine  and  its  derivatives,  MMH  (monomethylhydrazine) 
and  UDM  (unsymmetrical  dimethylhydrazine)  are  capable  of  exothermic 
decomposition,  and  can  be  used  as  monopropellants  as  well  as  the  fuel 
component  of  a  bipropellant  system.  In  either  case,  fuel  decomposition 
occurs  in  the  gas  phase  near  the  liquid  surface,  which  can  greatly  in¬ 
crease  the  drop  gasification  rate  in  comparison  to  conventional  bi¬ 
propellants.  When  burned  in  an  oxidizing  atmosphere,  the  decomposition 
flame  near  the  surface  can  be  surrounded  by  an  oxidization  zone,  yielding 
a  hybrid  structure.  Existing  information  on  the  combustion  mechanism 
of  hydrazine  fuels  was  relatively  limited,  in  view  of  the  complexities 
of  the  decomposition  process  and  the  possibility  of  hybrid  flame 
structures  (59-63) .  Therefore,  some  of  the  experimental  techniques 
developed  during  the  high  pressure  bipropellant  drop  combustion  study 
were  applied  to  the  hydrazine  fuels. 

A  second  area  of  interest  for  hydrazine  fuels  involves  the  problem 
of  combustion  instability.  Gasification  rates  for  monopropellant  drops 
vary  more  rapidly  with  pressure  than  is  the  case  for  conventional  bi- 
propellants.  This  suggests  a  greater  tendency  for  combustion  instability 
with  these  fuels  (see  Ref.  64  and  references  cited  therein) . 

Both  the  problem  of  steady  combustion  and  the  combustion  response 
of  drops  and  sprays  of  hydrazine  fuels  were  considered  during  the 
investigation. 

3.2  Steady  Drop  Combustion 

Introduction.  Earlier  work  on  hydrazine  combustion  (59-63) 
was  extended  by  burning  hydrazine,  MMH  and  UDMH  droplets  in  a  com¬ 
bustion  gas  environment.  The  test  environment  included  ambient  oxygen 
concentrations  in  the  range  0-42%  by  mass,  and  ambient  gas  temperatures 
of  1660-2530  K,  at  atmospheric  pressure.  Both  supported  drops  and 
porous  spheres  were  employed  in  the  tests,  in  order  to  obtain  a  rela¬ 
tively  large  drop  diameter  range,  1.1-19.1  nan. 

A  simplified  theoretical  model  of  hybrid  combustion  was  developed, 
employing  an  approach  suggested  by  Williams  (64) . 

Further  details  concerning  this  portion  of  the  investigation  can 
be  found  in  Refs.  7,  8,  26,  27  and  42. 
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Apparatus.  The  experiments  employed  the  flat  flame  burner  appa¬ 
ratus,  illustrated  in  Fig.  6.  For  suspended  droplet  tests,  the  drops 
were  mounted  on  a  100  pm  diameter  quartz  fiber.  The  arrangement  of  the 
porous  sphere  tests  was  similar  to  that  pictured  in  Fig.  7. 

Shadowgraph  motion  pictures  were  analyzed  to  provide  the  drop 
diameter  variation  as  a  function  of  time  for  the  suspended  drop  tests. 
This  information  allows  the  drop  burning  rate  to  be  computed  (26). 

The  fuel  flow  rate  to  the  sphere,  with  the  sphere  fully  wet  and  not 
dripping,  provides  the  burning  rate  directly,  for  the  porous  sphere 
tests. 

Theory.  The  objective  of  the  theoretical  model  was  to  provide  a 
semi-empirical  correlation  of  the  burning  rate  data.  The  model  con¬ 
siders  the  presence  of  both  a  decomposition  and  an  oxidation  flame. 

Each  reaction  zone  is  assumed  to  be  infinitely  thin.  The  position  of 
the  decomposition  flame  is  fixed  by  its  laminar  flame  speed,  which  is 
correlated  empirically.  The  position  of  the  oxidation  flame  is  fixed 
by  diffusion  rates  and  the  requirements  of  stoichiometry.  Convection 
effects  are  considered  using  a  film  theory  approximation.  The  same 
empirical  correlation  for  convection  used  for  the  bipropellant  com¬ 
bustion  theory  was  employed  to  compute  the  film  thickness.  The  low 
pressure  phase  equilibrium  model  is  adequate  for  the  test  conditions 
and  was  employed  in  the  analysis. 

The  theory  provides  predictions  of  the  position  of  the  two  flames, 
the  droplet  gasification  rates,  the  temperature  of  the  liquid  surface 
and  both  flame  temperatures . 

Results  and  Discussion.  The  burning  rate  of  hydrazine  as  a  function 
of  drop  diameter  is  illustrated  in  Fig.  12  for  a  typical  oxidation 
condition.  Burning  rate  estimations  from  both  the  hybrid  combustion 
model  and  the  conventional  bipropellant  model,  which  neglects  the 
presence  of  the  decomposition  flame,  are  also  shown  on  the  figure.  The 
bipropellant  results  are  shown  as  a  band  rather  than  a  single  curve,  in 
order  to  illustrate  the  effect  of  property  uncertainties.  The  hybrid 
model  contains  the  activation  energy,  E,  as  a  parameter  in  the  laminar 
flame  speed  correlation.  This  parameter  has  little  influence  on  the 
results  over  most  of  the  data  range  illustrated  in  Fig.  12.  The 
hybrid  model  clearly  provides  better  agreement  with  the  measurements. 
The  effect  of  reaction  is  illustrated  by  the  relatively  rapid  increase 
in  drop  burning  rate  with  size.  At  smaller  sizes,  however,  the  decom¬ 
position  flame  tends  to  merge  with  the  oxidation  flame,  and  the  two 
predictions  are  essentially  the  same.  Similar  results  were  obtained 
for  various  ambient  oxygen  concentrations,  temperatures  and  fuels  (25). 

Figure  13  is  an  illustration  of  the  comparison  between  all  the 
predictions  and  measurements  for  hydrazine.  In  addition  to  the  results 
of  the  present  investigation,  data  from  Refs.  59-62  are  also  considered. 
The  average  error  between  predictions  and  measurements  is  20  percent. 
Comparisons  for  MH  and  UDNH  yielded  similar  average  errors.  The  data 
range  for  the  comparison  is  summ  ed  in  Table  3. 


MASS  BURNING  RATE  (gm/sec) 
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DROP  DIAMETER  (cm) 

Fig.  12  Hydrazine  burning  rate  as  a  function  of  drop  diameter 
for  a  typical  oxidation  condition. 


EXPERIMENTAL  MASS  BURNING  RATE  (gm/sec) 
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PREDICTED  MASS  BURNING  RATE  (gm/sec) 

Fig.  13  Comparison  of  predicted  and  measured  burning  rates  for 
hydrazine. 


Table  3 


Summary  of  Test  Conditions  for  Drop  Combustion  Rates 
of  Hydrazine  Fuels 

Fuels  Hydrazine,  NWH,  UDMH 

Drop  diameters  0.38  -  19.1  mm 

Ambient  temperature  300-2530  K 

Ambient  oxygen  concentration  0-100% 


Pressure 


1-7.8  atm 
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Conclusions .  Although  a  more  rigorous  theory  would  be  desirable, 
the  comparison  between  the  hybrid  combustion  theory  and  the  experi¬ 
ments  was  satisfactory,  with  average  errors  less  than  20%.  Hydrazine, 

UDMH  and  MMH  all  exhibited  hybrid  combustion  effects  at  atmospheric 
pressure,  although  the  effect  was  most  pronounced  for  hydrazine.  The 
influence  of  pressure  has  not  been  clearly  resolved  during  this  investi¬ 
gation,  although  a  second-order  dependence  for  the  laminar  flame  speed 
of  hydrazine  is  suggested  by  the  measurements  of  Kosvic  and  Breen  (62) , 
at  least  at  pressures  up  to  7.8  atm.  Results  discussed  in  the  next 
section  confirm  the  second-order  dependence  at  pressures  greater  than 
atmospheric  pressure. 

3.3  Combustion  Response  of  Liquid  Hydrazine  Fuels 
3.3.1  Liquid  Strand  Combustion 

Introduction .  In  response  to  combustion  instability 
problems  in  liquid  fueled  rocket  engines,  several  studies  have  been 
completed  concerning  the  open-loop  combustion  response  of  bipropellant 
droplets  (63) .  The  present  investigation  considered  this  problem  for 
monopropellants . 

The  investigation  was  limited  to  the  open-loop  response  to  pressure 
fluctuations,  since  earlier  work  had  shown  that  response  to  velocity 
fluctuations  was  small  (63) .  The  regime  that  was  treated  does  not  in¬ 
volve  the  low  frequency  range  where  both  liquid-  and  gas-phase  processes 
are  quasisteady.  The  results  are  limited  to  conditions  where  transient 
liquid  and  gas-phase  phenomena,  or  thermal  waves,  aTe  interacting  with 
the  combustion  process.  The  limiting  case  of  a  constant  area  liquid 
strand  was  considered  in  the  first  phase  of  the  study.  This  portion  of 
the  study  is  discussed  in  this  section.  Strand  combustion  is  repre¬ 
sentative  of  very  large  droplets,  with  the  reaction  zone  close  to  the 
surface,  where  curvature  effects  are  small.  The  response  of  smaller 
drop  sizes  and  sprays  is  discussed  in  Section  3.3.2. 

The  predictions  were  verified  directly  by  experiment  -  apparently 
the  only  case  where  this  has  been  done  for  a  liquid  combustion  response 
problem.  Hydrazine  was  used  as  the  test  fuel,  since  it  can  sustain 
combustion  at  very  low  burning  rates.  With  such  a  material,  inter¬ 
actions  between  the  combustion  process  and  the  transient  liquid  thermal 
wave,  due  to  pressure  fluctuations,  can  be  observed  at  modest  frequencies, 
simplifying  the  measurements. 

Complete  details  concerning  this  portion  of  the  investigation  can 
be  found  in  Refs.  10,  12,  29  and  43. 

Apparatus .  A  sketch  of  the  apparatus  used  to  measure  the  open- loop 
response  of  liquid  hydrazine  strands  appears  in  Fig.  14.  The  liquid 
was  contained  in  a  glass  tube  and  ignited  at  the  top  with  an  electrically 
heated  wire.  Combustion  caused  the  liquid  surface  to  move  down  the 
tube,  past  the  observation  windows. 


Sketch  of  the  oscillatory  strand  combustion  apparatus 
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A  rotary  valve  arrangement  was  used  to  provide  oscillatory  pressure 
variations  within  the  test  chamber.  This  system  allowed  independent 
variation  of  the  mean  pressure,  and  the  amplitude  and  frequency  of  the 
pressure  oscillations. 

The  motion  of  the  liquid  surface  was  recorded  as  a  shadowgraph  on 
a  35  mm  strip  camera.  The  overall  optical  gain  of  the  surface  position 
measurements  was  12S.1. 

Theory.  The  analysis  for  the  open- loop  response  of  a  liquid  mono¬ 
propellant  was  similar  to  the  solid-propellant  combustion  analysis  of 
Tien  (65)  .  The  model  considers  both  gas-  and  liquid-phase  transient 
effects,  and  variable  properties  are  treated  in  the  gas  phase.  The 
combustion  process  is  assumed  to  be  one -dimensional,  premixed,  and 
laminar,  following  a  one-step,  irreversible  chemical  reaction.  In  the 
gas  phase,  the  molecular  weights  of  all  species  are  taken  to  be  constant 

and  equal,  the  gas  is  ideal,  the  Lewis  number  is  unity,  and  the  thermal 

conductivity  is  assumed  to  be  independent  of  composition  and  proportional 
to  temperature.  The  liquid  phase  properties  are  assumed  to  be  constant, 
with  negligible  reaction  and  dissolved  gas  concentrations.  Low  pressure 
phase  equilibrium  is  assumed  at  the  liquid  surface. 

A  perturbation  analysis  is  used  to  complete  the  solution.  Terms 
up  to  first-order  were  considered.  The  zero-order  results  were  checked 
with  steady  combustion  results,  prior  to  treating  combustion  response 
with  the  first-order  theory. 

Results  and  Discussion.  The  predicted  and  measured  steady  burning 
rates  of  hydrazine  as  a  function  of  pressure  are  illustrated  in  Fig.  15. 

Data  from  both  the  present  study  and  Antoine  (66)  are  shown.  The  two 

sets  of  measurements  are  in  reasonably  good  agreement.  The  present 
measurements  are  properly  corrected  for  flame  curvature,  due  to  surface 
tension,  and  probably  are  more  accurate  than  the  results  of  Ref.  66. 

The  theoretical  results  illustrated  in  Fig.  15  were  obtained  by 
fitting  the  zero-order  (steady)  theory  at  a  particular  condition,  for 
various  reaction  orders  and  activation  energies.  Subsequent  results 
are  predictions.  The  effect  of  activation  energy  on  the  correlation 
between  predictions  and  measurements  is  small.  The  reaction,  however, 
shifts  from  first  to  second  order  at  pressures  near  2  atm  and  two 
separate  correlations  are  needed  to  cover  the  range  of  the  data.  Pre¬ 
dictions  of  liquid  surface  temperatures  and  the  temperature  distribution 
in  the  liquid  phase  were  also  examined.  Agreement  between  theory  and 
experiment  was  excellent  (31,43). 

Having  obtained  a  reasonably  satisfactory  model  for  steady  com¬ 
bustion,  oscillatory  combustion  was  considered.  Since  the  theory  is 
linear,  the  effect  of  the  amplitude  of  the  pressure  oscillation  was 
examined.  It  was  found  that  the  amplitude  of  the  burning  rate  oscil¬ 
lation  was  linear  for  pressure  oscillations  less  than  20  percent  of 
the  mean  pressure.  Subsequent  tests  were  limited  to  pressure  ampli¬ 
tudes  less  than  15  percent  of  the  mean  pressure,  in  order  to  assure 
operation  in  the  linear  range  for  comparison  with  the  linear  response 
theory. 
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Fig.  15  Theoretical  and  experimental  steady  hydrazine  strand 
burning  rates  as  a  function  of  pressure. 
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Figures  16  and  17  are  illustrations  of  the  comparison  between 
predicted  and  measured  amplitudes  and  phase  angles  of  the  liquid 
surface  oscillations  as  a  result  of  imposed  pressure  oscillations  of 
a  given  frequency.  Theoretical  results  are  illustrated  for  various 
activation  energies.  A  second-order  reaction  is  assumed  for  the  mean 
pressure  of  these  results,  9.77  atm,  c.f.  Fig.  15.  The  best  match 
between  theory  and  experiment  is  obtained  for  a  dimensionless  acti¬ 
vation  energy  in  the  range  10-15. 

Employing  the  chemical  parameters  determined  from  the  experiments, 
the  real  part  of  the  acoustic  admittance  was  computed  for  hydrazine 
strand  combustion.  A  typical  result  is  illustrated  in  Fig.  18.  As  the 
imposed  frequency  is  increased,  two  positive  admittance  peaks  are  ob¬ 
served.  The  first  peak  is  due  to  the  interaction  between  the  combustion 
process  and  liquid  phase  transient  effects  while  the  second  peak  is  due 
to  gas  phase  transient  processes.  These  characteristics  are  similar 
to  T'ien's  (65)  results  for  solid  propellants.  Increasing  the  pressure 
tends  to  increase  the  magnitude  of  the  first  peak  and  decrease  the 
magnitude  of  the  second  peak. 

Combustion  instability  is  usually  associated  with  values  of  the 
real  part  of  the  acoustic  admittance  greater  than  a  value  on  the  order 
of  unity.  At  atmospheric  pressure,  only  the  gas  phase  transient  range 
has  sufficient  amplification  to  drive  instability.  At  pressures  greater 
than  10  atm,  however,  both  regimes  are  significant.  Table  4  is  a  sum¬ 
mary  of  the  actual  frequency  ranges  where  the  acoustic  admittance  exceeds 
unity.  Combustion  instability  is  usually  associated  with  frequencies 
in  the  range  500-30,000  Hz  (63).  At  pressures  greater  than  10  atm, 
the  ranges  shown  in  Table  4  suggest  that  the  liquid  phase  mode  is 
the  dominant  regime  in  this  frequency  range. 

Conclusions.  The  combustion  mechanism  of  hydrazine  shifts  from 
first  to  second-order  at  pressures  slightly  greater  than  atmospheric 
pressure.  Linear  combustion  response  was  observed  directly,  the  only 
case  where  this  has  been  done  for  liquids.  Theory  and  experiment  were 
in  good  agreement  with  respect  to  both  mean  properties  and  combustion 
response.  The  analysis  correctly  predicted  a  complex  variation  of  the 
phase  shift  with  pressure.  Both  liquid  and  gas  phase  transient  effects 
interact  with  the  combustion  process  to  produce  frequency  ranges  when 
the  acoustic  admittance  exceeds  unity.  However,  liquid  transient 
phenomena  appear  to  be  more  important  for  the  stability  of  typical 
combustion  systems  at  pressures  greater  than  10  atm.  These  results 
are  limited  to  burning  liquid  hydrazine  strands,  which  are  only  repre¬ 
sentative  of  the  largest  drops  in  actual  combustion  chambers. 

3.3.2  Drop  and  Spray  Combustion 

Introduction.  Large  monopropellant  drops  behave  similar  to 
burning  liquid  strands,  and  would  show  positive  combustion  response 
regions  similar  to  the  results  discussed  in  Section  3.3.1.  On  the 
other  hand,  reaction  effects  become  small  for  small  drops,  suggesting 


AMPLITUDE  OF  SURFACE  OSCILLA 


Fig.  16  Theoretical  and  experimental  magnitude  of  the  liquid 

surface  oscillation  of  burning  hydrazine  as  a  function 
of  the  frequency  of  the  imposed  pressure  oscillation. 
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FREQUENCY  (Hz) 


Fig.  17  Theoretical  and  experimental  phase  angle  of  the  liquid 
surface  oscillation  of  burning  hydrazine  as  a  function 
of  the  frequency  of  the  imposed  pressure  oscillation. 


Fig,  18  Predicted  real -part  of  the  acoustic  admittance  for  hydrazine  strand 
combustion. 


Table  4 


Frequency  Ranges  for  Potentially  Unstable 
Combustion  of  Hydrazine 


Mean  Pressure 
(atm) 

Liquid  Transient 
Range  (Hz) 

Gas  Transient 
Range  (Hz) 

1 

- 

300-1500 

10 

15  (approx) 

4500-21500 

100 

450-1650 

>225,000 

37 


that  their  response  characteristics  would  be  similar  to  conventional 
bipropellant  drops.  Earlier  results  indicated  that  the  response  of 
bipropellant  drops ,  by  interaction  of  combustion  with  liquid  and  gas 
phase  transient  effects,  was  small  (63).  The  objective  of  the  present 
investigation  was  to  examine  the  transition  between  these  two  regimes 
and  to  also  make  a  preliminary  assessment  of  spray  response. 

The  study  was  limited  to  hydrazine  combustion.  The  theory  was 
largely  based  on  the  earlier  results  for  hydrazine  strand  burning 
response.  As  noted  earlier,  the  response  under  consideration  is  not 
that  of  a  combustion  chamber  as  a  whole,  where  both  liquid  and  gas 
phases  are  quasisteady.  Emphasis  was  placed  on  the  interaction  between 
the  combustion  process  and  the  transient  thermal  wave  in  the  liquid 
phase,  in  the  presence  of  imposed  pressure  oscillations. 

Further  details  on  this  portion  of  the  study  may  be  found  in 
Refs.  14,  15,  32,  33  and  45. 

Theory.  Most  of  the  basic  elements  of  the  model  were  the  same  as 
for  the  strand  combustion  case.  The  model  considers  gas-and  liquid- 
transient  effects.  Variable  physical  properties  are  treated  in  the 
gas  phase.  To  simplify  the  analysis,  forced  convection  was  neglected, 
providing  a  spherically  symmetric  combustion  process.  This  limitation 
is  less  serious  for  monopropellants  than  bipropellants.  The  response 
is  greatest  where  the  reaction  zone  is  near  the  liquid  surface,  a 
condition  where  the  combustion  process  is  not  strongly  affected  by 
forced  convection  (67,68). 

The  ambient  pressure  is  assumed  to  be  oscillating  with  a  wavelength 
that  is  long  in  comparison  to  the  dimensions  of  the  drop  combustion 
field.  The  period  of  oscillation,  however,  is  assumed  to  be  short  in 
comparison  to  the  drop  lifetime.  This  assumption  applies  a  low 
frequency  bound  for  the  analysis.  Low  pressure  phase  equilibrium  was 
assumed  at  the  liquid  surface. 

A  common  assumption  for  drop  combustion  studies  is  that  the  drop 
is  heated  uniformly  to  its  wet-bulb  temperature.  It  is  generally 
recognized,  however,  that  temperature  gradients  persist  within  drops 
throughout  their  lifetime,  as  artifacts  of  the  heat-up  period.  The 
high  mean  combustion  rate  of  monopropellants  increases  the  possibility 
of  temperature  gradients  (69) .  In  order  to  investigate  the  potential 
affect  of  temperature  gradients  on  combustion  response,  two  limiting 
cases  were  treated:  i)  mean  liquid  temperature  varying  as  though  a 
quasisteady  combustion  wave  was  propagating  toward  the  center  of  an 
unheated  drop,  and  ii)  mean  liquid  temperature  constant  at  the  wet- 
bulb  temperature. 

Similar  to  the  strand  analysis,  a  perturbation  solution  was  em¬ 
ployed.  The  zero-order  solution  provides  the  steady  combustion  charac¬ 
teristics.  The  first-order  solution  provides  the  linear  combustion 
response . 
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Analysis  was  also  completed  to  determine  the  response  of  a  spray, 
given  the  response  of  individual  drops.  This  response  was  defined  in 
terms  of  the  mass  burning  rate  fluctuation,  averaged  over  the  particular 
drop  size  distribution  in  the  spray,  normalized  by  the  total  steady 
burning  rate  of  the  spray. 

Results  and  Discussion.  Figure  19  is  an  illustration  of  droplet 
response  as  a  function  of  frequency  and  drop  size  (the  parameter  A  is 
proportional  to  drop  diameter  squared),  for  non-wet -bulb  combustion, 
l.ie  gas  phase  is  assumed  to  be  quasisteady  for  these  results.  The 
response  generally  increases  as  A,  or  the  drop  size,  increases.  For 
non-wet-bulb  combustion,  a  response  peak  is  observed  at  frequencies 
near  the  characteristic  frequency  for  liquid  phase  transient  effects 
(this  is  not  observed  for  wet-bulb  combustion).  For  large  values  of  A, 
the  system  approaches  the  characteristics  of  strand  combustion.  In  this 
range,  the  response  is  unity  at  low  frequencies,  where  both  phases  are 
quasisteady.  This  occurs  since  the  burning  rate  is  proportional  to 
pressure  for  steady  strand  combustion  with  a  second-order  reaction. 
Droplets  evaporating  with  little  reaction  near  the  surface  (low  values 
of  A)  generally  exhibit  low  or  negative  response  for  all  frequencies, 
in  agreement  with  earlier  studies  for  bipropellant  drop  response. 

Figure  20  is  an  illustration  of  spray  response  calculations  for 
non-wet-bulb  combustion.  Several  spray  distributions  are  considered 
involving  initial  injection  of  varying  proportions  of  different  drop 
size  groups  (ranging  from  100  percent  large  drops  to  100  percent  small 
drops) .  The  total  response  approaches  unity  when  che  two  largest 
sized  drops  comprise  only  about  10  percent  of  the  injected  mass, 
corresponding  to  0.3  percent  of  the  drop  injection  rate.  This  indi¬ 
cates  that  even  a  very  small  percentage  of  large  drops  in  a  spray  can 
result  in  substantial  positive  response.  This  occurs  since  the  large 
drops  comprise  an  appreciable  fraction  of  the  drops  within  a  combustion 
chamber,  due  to  their  long  lifetime,  e.g.,  when  the  two  largest  groups 
comprise  only  11  percent  of  the  injected  spray  mass,  drops  originating 
from  these  groups  provide  90  percent  of  the  mass  of  fuel  within  the 
burning  spray. 

Conclusions.  The  response  of  monopropellant  drops  and  sprays  can 
exceed  unity  for  drop  sizes  and  frequencies  representative  of  high- 
frequency  combustion  instability,  providing  a  viable  mechanism  for  this 
phenomena.  In  contrast,  drop  response  for  bipropellants  yields  low 
response  values.  The  response  is  greatest  for  large  drops,  having 
mean  liquid  temperature  gradients,  at  high  pressures.  Response  factors 
were  less  than  unity  for  drops  at  their  wet-bulb  state.  The  response 
of  a  monodisperse  spray  is  less  than  the  response  of  the  injected 
drop-size,  due  to  the  reduced  response  of  small  drop?.  Response  de¬ 
pends  very  strongly  on  the  initial  drop  size  distribution  for  poly- 
disperse  sprays;  however,  even  a  small  percentage  of  large  drops  can 
provide  appreciable  response  due  to  their  long  lifetimes  and  large 
response . 


Fig.  20  Predicted  response  of  a  polydisperse  hydrazine 
spray  as  a  function  of  frequency. 
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IV.  COMBUSTION  OF  SPRAYS 

4.1  Introduction.  The  development  of  turbulence  models  with 
second-order  closure,  and  the  availability  of  general-purpose  computer 
programs  to  solve  complex  systems  of  flow  and  combustion  equations, 
has  provided  the  potential  to  develop  spray  combustion  models  to 
assist  combustor  development  efforts.  In  fact,  models  along  these 
lines  are  already  being  used  by  some  gas  turbine  manufacturere  (36) . 

A  difficulty  encountered  with  spray  models,  however,  is  the  absence 
of  data  concerning  spray  characteristics  produced  by  a  given  injector 
(injected  drop  size  and  velocity  distributions  for  hot-firing 
conditions) . 

A  long  term  objective  of  spray  combustion  studies  is  to  develop 
comprehensive  models  where  all  factors  can  be  predicted,  including 
injector  characteristics.  There  are  many  difficulties  to  be  resolved, 
however,  and  simplified  models  can  still  provide  useful  results  in  the 
interim.  One  simplification  that  can  be  made,  involves  the  assumption 
of  locally  homogeneous  flow  (LHF) .  This  approximation  implies  that 
there  is  no  velocity  difference  between  the  phases  (no  slip)  and  that 
the  phases  are  in  thermodynamic  equilibrium,  i.e.,  that  the  temperature 
and  fugacity  of  all  components  is  the  same  in  all  phases.  This  approxi¬ 
mation  is  most  appropriate  for  finely  atomized  sprays,  where  individual 
drops  can  rapidly  respond  to  changes  in  their  local  environment,  i.e., 
cases  where  the  process  is  mixing  controlled  as  opposed  to  being 
vaporization  controlled.  All  real  sprays  exhibit  some  degree  of  slip 
and  loss  of  thermodynamic  equilibrium.  In  this  case,  the  results  of 
LHF  ,  if  accurate,  represent  an  upper  bound  to  the  rate  of  development 
of  the  combustion  process  -  the  actual  spray  will  develop  more  slowly. 

Locally  homogeneous  flow  models  are  widely  used  in  the  analysis 
of  two-phase  heat  transfer  and  hydraulic  problems  (70) .  Models  of 
this  type  have  also  been  employed  for  gas-liquid  jets  with  success, 
since  bubbles  respond  rapidly  to  surrounding  liquid  flows  (71,72). 

There  has  also  been  earlier  application  of  the  method  to  sprays.  Onuma 
and  Ogasawara  (73)  demonstrate  startling  similarities  between  measured 
properties  in  gas  flames  and  well  atomized  sprays.  Newman  and  Bruzustowski 
(74)  also  successfully  apply  the  concept  to  sprays  evaporating  near  their 
thermodynamic  critical  point.  Models  along  these  lines  have  also  been 
successfully  applied  to  Diesel  combustion  (36).  On  the  other  hand, 

Khalil  and  Whitelaw  (75)  evaluated  a  gas  simulant  model  (this  is  not 
a  complete  locally  homogeneous  flow  model,  only  the  energy  release, 
momentum  and  stoichiometry  is  matched,  not  the  complete  equation  of 
state  as  in  traditional  locally  homogeneous  flow  models)  with  less 
success,  finding  significant  discrepancies  between  predictions  and 
measurements  for  sprays  having  a  Sauter  mean  diameter  as  small  as 
45  pm.  Thus,  adequate  criteria  for  the  validity  of  LHF  models  had  not 
been  determined  in  earlier  work. 
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The  present  investigation  extends  this  work,  attempting  to  provide 
a  more  systematic  evaluation  of  the  LHF  model.  A  systematic  method  of 
estimating  the  conditions  where  LHF  models  are  accurate  was  also 
sought.  The  study  was  limited  to  a  finely  atomized  spray,  injected 
into  a  stagnant  environment,  in  order  to  simplify  the  flow  field.  The 
theory  employs  a  k-e-g  turbulence  model,  similar  to  those  used  for 
gaseous  diffusion  flames  (76) .  This  model  has  the  advantage  of 
minimizing  the  number  of  empirical  parameters  to  be  specified. 

Since  turbulence  models  contain  a  large  number  of  empirical 
parameters,  the  evaluation  proceeded  through  a  series  of  simpler 
flows  before  sprays  were  considered.  The  sequence  was  as  follows: 

i.  constant  density  jets 

ii.  variable  density  jets 

iii .  evaporating  sprays 

iv.  combusting  jets 

v.  combusting  sprays 

The  entire  investigation  was  divided  into  two  phases.  The  first  phase 
considered  flows  i-iii,  while  the  second  phase  considered  iv  and  v. 

Each  of  these  phases  is  discussed,  in  turn,  in  the  following. 

Further  details  concerning  this  portion  of  the  investigation  can 
br  found  in  Refs.  16-18,  35-38,  and  46-49. 

4.2  Evaporating  Sprays 

Introduction.  In  this  portion  of  the  study,  an  experimental 
apparatus  was  constructed  in  order  to  conduct  measurements  within 
isothermal  air  jets,  dense  gas  jets  and  an  evaporating  spray.  A 
variety  of  measurements  were  completed  for  comparison  with  predictions. 
Jet  flow  rates  and  the  total  momentum  of  the  jet  were  measured  to 
provide  initial  conditions.  Profiles  of  mean  velocity,  velocity 
fluctuations  and  Reynolds  stress  were  measured  for  all  the  flows.  Mean 
concentration  measurements  were  completed  for  the  isothermal  gas  jet 
and  the  evaporating  spray.  Mean  temperature  and  dropsize  distribution 
measurements  were  completed  for  the  evaporating  spray.  Both  theoretical 
and  experimental  methods  used  in  this  study  were  checked  against  earlier 
work  in  the  literature  (76-81) . 

The  major  elements  of  the  theory  were  also  developed  in  this  phase 
of  the  study.  Due  to  the  approach  taken  in  the  analysis,  each  specific 
case  is  only  distinguished  by  a  new  equation  of  state.  Thus  subsequent 
considerations  of  combusting  flows  involve  little  additional  analysis. 

Apparatus .  All  testing  was  done  in  the  same  apparatus.  The  in¬ 
jector  (air  atomizing  nozzle.  Spraying  Systems  Company  Model  1/4  J  2050 
with  a  No.  67147  air  nozzle)  was  oriented  vertically  downward.  The 
test  area  was  2  m  high  by  0.92  m  square,  enclosed  with  16  mesh  screen. 
Two  linear  transversing  mechanisms  were  used  to  position  the  nozzle  in 
a  horizontal  plane.  The  nozzle  could  also  be  moved  in  the  vertical 
direction,  along  slides. 
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Gas  flow  rates  through  the  injector  were  metered  with  critical  flow 
orifices.  Freon  11  was  used  as  the  spray  liquid,  its  flow  rate  was 
measured  directly  by  observing  the  liquid  level  displacement  in  the 
supply  vessel.  The  liquid  flow  rate  was  varied  by  changing  the  air 
pressure  over  the  surface  of  the  liquid.  The  momentum  of  the  jet  was 
determined  directly  by  impingement  on  a  flat  plate,  supported  by  a 
load  cell. 

Velocities  and  turbulence  quantities  were  measured  by  a  single 
channel  laser-Doppler  anemometer.  The  arrangement  of  the  LDA  is  illus¬ 
trated  in  Fig.  21.  The  single  channel  instrument  was  employed  for 
various  velocity  components  and  turbulence  quantities  by  changing  the 
orientation  of  the  beams  (17,47).  Temperatures  were  measured  with  a 
25  pm  00  wires  mounted  in  crossflow  from  heavier  leads  spaced  20  mm 
apart.  Overall  compositions  were  obtained  by  isokinetically  sampling 
the  flow  with  a  heated  probe  and  analyzing  the  sample  using  a  gas 
chromatograph.  Liquid  mass  flux  and  the  drop  size  distribution  in 
sprays  were  obtained  using  a  particle  impaction  technique. 

Three  test  conditions  were  considered,  an  air  jet  (constant 
density  jet),  an  SF6  jet  (variable  density  jet)  and  an  evaporating 
Freon  11  spray.  The  test  conditions  are  summarized  in  Table  5. 

Theory.  The  theoretical  model  assumes  an  axisymmetric  turbulent 
jet  in  an  infinite  stagnant  gaseous  environment.  The  spray  is  treated 
with  LHF  approximation.  This  implies  that  drop  and  gas  velocities  and 
temperatures  are  equal  at  each  point  in  the  flow,  and  that  the  drops 
and  the  vapor  are  in  thermodynamic  equilibrium,  i.e.,  that  the  partial 
pressure  of  the  vapor  is  the  saturated  vapor  pressure  of  the  liquid, 
as  long  as  any  liquid  is  present.  The  calculations  are  based  on  Reynolds 
averaging  the  conservation  equations.  Similar  to  other  models  of  this 
type  (76),  production  and  dissipation  effects  due  to  density  fluctuations 
are  ignored.  It  is  also  assumed  that  laminar  transport  occurs  with 
equal  diffusion  coefficients,  that  turbulent  species  and  heat  exchange 
coefficients  are  equal  as  well,  and  that  radiation  and  kinetic  heating 
are  negligible. 

Under  these  assumptions,  the  conservation  equations  are  as  follows 
(17,76): 


7 1? 


D  (?)  *  0 


(1) 


D  (u)  -  a  (p^-p) 


(2) 


D  (k) 


(3) 
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Table  5 

Summary  of  Test  Conditions3 


Case 

1 

2 

3 

Injector  Fluid 

Air 

Sulfur  Hexafluoride 

Freon- 11 

Injector  Flow  Rates 
(g/s) 

Gas 

0.467 

0.103 

0.225 

Liquid 

0 

0 

1.548 

Injector  Pressures 
(kPa) 

Gas 

203 

186 

211 

Liquid 

— 

— 

138 

Jet  Momentum 

(N) 

0.100 

0.0093 

0.132 

Sauter  Mean  Diameter*3 
(pm) 

— 

— 

29 

Initial  Velocity 
(m/s) 

214.1 

90.30 

74.45 

3T00  -  296  K,  Pa  -  97  kPa 
^Measured  203  mm  from  injector. 
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where 


D  Cc)  '  cel  “t  I  (3?>  *  ce2S  f“  <4> 


D  <*>  ■  Vt  -  cg25  F  (5) 


p  v°  =  p  v  +  p'v'  (6) 


D(4»  -  P  G  ||  ♦  p  v°  H  -  I  ^  (r  ^  ||)  (7) 


for  <p  *  u,  i,  k,  e  or  g.  The  boundary  conditions  for  these  equations 
are: 


r  -  0,  If  =  0;  r  -  «,  *  =  o 


(8) 


The  initial  conditions  are: 


x  »  0,  r  <  d/2,  u  »  M0/m0>  f  a  1,  g  *  0,  k  »  kQ,  e  *  eQ  (9) 

where  the  parameters  kQ  and  eQ  are  assumed  to  have  values  appropriate 
for  fully  developed  flow  in  tubes.  The  turbulent  viscosity  is  obtained 
from  the  following  constitutive  equation: 


ut  *  cp  p  *2/e 


(10) 


The  model  is  completed  by  relating  the  mixture  fraction  to  other 
properties  of  the  flow,  p,  Y. ,  T,  etc.  Since  f  is  a  passive  scalar, 
local  properties  are  completely  determined  by  the  thermodynamic 
equilibrium  state  attained  during  adiabatic  mixing  of  an  amount  f  of 
injector  fluid  and  (1-f)  of  ambient  fluid,  at  the  ambient  pressure. 
Determination  of  the  equation  of  state  then  involves  solution  of 
equations  for  conservation  of  mass  and  energy  for  this  mixture  along 
with  and  phase  equilibrium  requirements  if  the  mixture  is  not  single¬ 
phase. 


4>- 
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Figure  22  is  an  illustration  of  the  equation  of  state  for  the 
Freon  11  spray  considered  in  the  experiments.  In  this  case,  both  air 
and  Freon  11  leave  the  injector,  while  the  ambient  gas  is  pure  air. 

As  the  mixture  fraction  decreases  from  unity,  the  concentration  of 
Freon  liquid  decreases,  with  all  the  liquid  disappearing  at  a  mixture 
fraction  of  0.33.  At  this  point,  the  temperature  is  a  minimum  and  the 
concentration  of  Freon  11  vapor  is  a  maximum.  The  presence  of  liquid 
Freon  11  results  in  a  rather  nonlinear  density  variation,  particularly 
at  high  mixture  fractions. 

Given  the  equation  of  state,  the  mean  value  of  any  scalar  quantity 
(other  than  ?,  k,  c  and  g)  can  be  determined  if  the  probability  density 
function  for  f  is  known  as  a  function  of  position  in  the  flow.  The 
mean  value  is  then  given  by 


$  *  *(f)  P  (£)  df  (11) 

Jo 


where  p(f)  is  the  probability  density  function  for  f.  Following  Lock- 
wood  and  Naguib  (76),  a  clipped  Gaussian  probability  density  function 
was  assumed.  There  are  two  parameters  in  this  function  which  must  be 
specified,  the  most  probable  value  of  the  distribution  and  the  variance. 
These  parameters  are  found  as  a  function  of  position  by  noting 

?  »  [1fp(f)df,  g  *  f  (f-f)2p(f)df  (12) 

J0  J0 


Both  ?  and  g  are  known  from  the  integration  of  Eqs.  (1) - (10) ,  providing 
two  equations  to  solve  for  the  unknown  parameters  in  the  probability 
density  distribution. 

The  calculations  were  performed  using  a  modified  version  of 
general-purpose  parabolic  flow  solver,  GENMIX  (82).  The  enpirical 
constants  appearing  in  the  model  were  the  same  as  those  recommended 
by  Lockwood  and  Naguib  (76).  The  values  are  summarized  in  Table  6. 

These  constants  were  used  in  all  regions  of  the  flow,  including  the 
shear  layer  near  the  injector  exit. 

Results  and  Discussion.  Figure  23  is  an  illustration  of  the 
radial  variation  of  mean  axial  velocity  for  an  isothermal  single-phase  jet, 
a  variable  density  single-phase  jet  and  the  evaporating  spray.  The  axial 
velocity  is  normalized  by  the  centerline  velocity.  The  radial  distance 
is  normalized  by  the  axial  distance  from  the  injector  so  that  a  direct 
measure  of  the  width  of  the  flow  can  be  obtained.  Experimental  results 
from  Wygnanski  and  Fiedler  (78)  and  Hetsroni  and  Sokalov  (80)  are 


MASS  FRACTION  p(kg/m3)  T(K) 


Table  6 


Constants  In  the  Turbulence  Model 


Constant 

Value 

cw 

0.09 

c 

1.44 

el 

c  -  c 
e2  g2 

1.89,  1.84a 

C 

2.8 

gl 

ak 

1.0 

a 

1.3 

e 

af 

0.7 

CJ 

_S _ 

0.7 

g 

Constant  density  and  variable  density  flows, respectively. 
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shown,  along  with  the  present  measurements.  The  theoretical  results 
for  each  case  also  appear,  for  comparison  with  the  measurements.  In 
general,  the  comparison  between  the  measurements  of  the  various 
studies,  and  between  predictions  and  measurements,  is  excellent. 

Similar  agreement  was  also  obtained  for  mixture  fraction. 

Figure  24  is  an  illustration  of  the  radial  variation  of  Reynolds 
stress  for  the  three  flows.  The  Reynolds  stress  is  normalized  by  the 
square  of  the  centerline  velocity.  The  radial  coordinate  is  the  same 
as  Fig.  23.  The  measurements  are  from  Wygnanski  and  Fiedler  (78)  and 
the  present  investigation.  Theory  is  also  illustrated  for  comparison 
with  the  measurements.  The  scatter  of  data  is  greater  for  this  turbu¬ 
lence  quantity  than  for  mean  quantities.  However,  the  general  agree¬ 
ment  between  predictions  and  measurements  is  satisfactory. 

Greater  differences  between  the  flows  become  apparent  when  the 
axial  variation  of  mean  velocity  and  mixture  fraction  are  considered. 
Figures  25  and  26  are  illustrations  of  these  parameters,  as  a  function 
of  distance  from  the  injector.  In  addition  to  the  present  data,  re¬ 
sults  from  other  investigations  of  single  and  two-phase  jets  are 
shown  (77-79,  81).  Earlier  experimental  results  for  air  jets  are  in 
good  agreement  with  the  present  measurements.  Predictions  of  the 
theory  are  also  illustrated  for  each  case.  The  single-phase  jets  are 
modeled  quite  well,  including  variable  density  effects.  The  dis¬ 
crepancies  are  larger  for  the  two-phase  jets.  For  the  air  jet  into 
water,  the  velocity  variation  is  predicted  quite  well,  but  the  error 
is  significant  for  mixture  fraction.  The  method  used  by  Tross  (81)  to 
evaluate  mixture  fraction,  tends  to  underestimate  this  quantity  and  the 
problem  could  be  with  the  experiments.  This  flow  also  has  the  largest 
density  variation  of  all  the  flows  considered,  and  a  difficulty  with 
neglecting  density  fluctuations  in  the  model  cannot  be  ruled  out 
(17,37).  The  measured  variation  of  mean  velocity  and  mixture  fraction 
for  the  evaporating  spray  clearly  lags  the  predictions.  Comparisons 
between  predicted  and  measured  liquid  fluxes  and  temperatures  also  in¬ 
dicated  that  the  flow  was  developing  more  slowly  than  predicted. 

Since  the  other  flows  were  predicted  reasonably  well  by  the  model, 
it  appeared  that  slip  and  loss  of  thermodynamic  equilibrium  might  be 
the  reason  for  the  errors  in  the  spray  predictions.  This  was  examined 
by  completing  drop-life-histories  in  the  flow.  The  life-history  cal¬ 
culations  considered  drops  passing  along  the  centerline  of  the  flow 
with  initial  conditions  corresponding  to  properties  at  the  injector 
exit.  The  local  ambient  conditions  of  the  drop  were  obtained  from  the 
predictions  of  the  LHF  turbulence  model.  Since  computation  of  drop- 
life-histories  is  subject  to  errors  due  to  uncertainties  in  the 
properties  to  be  used  in  the  computations,  the  calculation  procedure 
was  calibrated  using  experimental  results  from  single  drop  experi¬ 
ments  with  well-defined  ambient  conditions  around  the  test  drop. 

The  results  of  the  drop-life-history  calculations  are  illustrated 
in  Fig.  27.  Initial  drop  diameters  of  10,  30  and  50  pm  are  considered. 
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The  diameter,  velocity  and  temperature  of  the  drop,  and  the  velocity 
and  temperature  of  the  gas  at  the  centerline  of  the  spray,  are  given  as 
a  function  of  position  from  the  injector.  The  drop  velocities  clearly 
decrease  less  rapidly  than  the  gas  velocities.  Appreciable  levels  of 
slip  occur  in  the  region  just  downstream  of  the  potential  core  of  the 
jet.  There  are  also  significant  differences  in  the  temperatures  of  the 
two  phases.  The  temperature  errors  are  largest  in  the  downstream 
region,  where  the  gas  temperature  tends  to  increase  while  the  drop 
temperature  tends  to  approach  a  wet-bulb  value.  The  life-history  cal¬ 
culations  indicate  that  drops  are  present  for  values  of  x/d  greater 
than  500,  which  agrees  with  experimental  observations.  It  appears 
that  the  spray  considered  in  this  investigation  would  have  to  have 
initial  drop  diameters  less  than  10  pm,  in  order  to  obtain,  quantitative 
accuracy  with  the  LHF  model. 

Conclusions.  The  model  provides  excellent  predictions  of  mean 
quantities  in  single-phase  constant  and  variable  density  jets.  The 
agreement  between  predicted  and  measured  turbulence  quantities  in 
single-phase  jets  was  satisfactory.  The  model  provided  useful 
qualitative  results  for  the  spray,  but  slip  and  loss  of  thermodynamic 
equilibrium  limited  the  quantitative  predictions  of  the  LHF  model, 
even  for  a  spray  having  a  Sauter  mean  diameter  of  29  pm  at  atmospheric 
pressure  (maximum  errors  were  in  the  range  20-40%  for  axial  variations 
of  mean  velocity  and  mixture  fraction) .  Analysis  of  drop-lift  histories 
suggests  that  a  Sauter  mean  diameter  less  than  10  pm  would  be  required 
for  quantitative  accuracy  at  these  conditions.  These  small  sizes  are 
needed  due  to  the  relatively  rapid  development  of  the  flow,  resulting 
from  the  relatively  small  injector  diameter. 

Since  injector  diameters  are  generally  small,  this  finding  implies 
that  accurate  analysis  of  most  practical  sprays  will  require  models 
which  allow  for  slip  and  finite  rate  processes  between  the  phases.  The 
LHF  model  is  still  of  value  in  cases  where  the  necessary  information  on 
spray  characteristics  is  limited,  reducing  the  potential  accuracy  of 
slip  models.  There  are  relatively  few  parameters  to  be  prescribed  for 
the  LHF  model,  other  than  a  conventional  thermodynamic  prediction  for 
the  equation  of  state.  The  results  provide  a  lower  bound  for  spray 
operation  -  since  actual  sprays  tend  to  develop  more  slowly  than  the 
LHF  predictions.  Conducting  drop  life-history  calculations,  using  the 
LHF  model  to  prescribe  the  drop  environment,  is  a  convenient  means  of 
evaluating  the  potential  accuracy  of  LHF  predictions. 

4.3  Spray  Combustion 

Introduction.  The  final  stage  of  this  investigation  con¬ 
sidered  the  application  of  the  LHF  model  to  combusting  sprays.  Two 
cases  were  considered,  a  spray  from  an  air-atomizing  injector  at 
atmospheric  pressure  and  a  spray  from  a  pressure  atomizing  injector  at 
high  pressures.  The  former  case  is  a  potential  application  of  the  LHF 
model,  due  to  the  small  drop  diameters  produced  with  air  atomizing 
injectors  and  the  relatively  large  flame  length  at  atmospheric  pressure 
(similar  to  the  experimental  conditions  of  Onuma  and  Ogasawara  (73)). 

LHF  models  also  show  promise  at  high  pressures,  since  the  densities  of 
the  two  phases  approach  one  another.  The  LHF  method  also  provides  a 
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relatively  simple  approach  for  treating  combustion  near  the  thermo¬ 
dynamic  critical  point,  where  difficulties  are  encountered  with 
drop- life-history  computations. 

Apparatus .  Tests  at  atmospheric  pressure  were  conducted  in  the 
apparatus  discussed  in  Section  4.2.  The  major  difference  was  that 
the  injector  faced  upward.  Test  conditions  involved  repeating  the 
air  and  sulfurhexafluoride  tests,  discussed  earlier;  and  conducting 
tests  in  a  combusting  propane  gas  jet  and  a  combusting  pentane-fueled 
spray.  The  flames  were  stabilized  at  the  injector  exit  using  four 
small  capillary  flames  fueled  with  hydrogen.  Measurements  in  the 
combusting  jets  were  limited  to  mean  axial  velocities,  using  the  LDA, 
and  gas  temperatures,  using  50  pm  diameter  Pt/Pt-10%  Rh  thermocouples. 
The  thermocouple  wires  were  shielded  in  the  spray  flame,  in  order  to 
prevent  impingement  of  drops  on  the  junction.  The  test  conditions 
for  the  combusting  jets  at  atmospheric  pressure  are  summarized  in 
Table  7. 

Tests  at  elevated  pressures  were  conducted  in  the  apparatus 
pictured  in  Fig.  28.  The  environment  of  the  spray  consisted  of  air, 
under  pressure,  within  the  test  chamber.  Since  only  a  limited 
quantity  of  air  was  present,  the  injector  was  operated  for  a  relatively 
short  time,  in  order  to  reduce  vitiation  of  its  atmosphere.  The  test 
duration  was  long  enough,  however,  to  insure  steady  combustion. 

Liquid  fuel  for  the  test  was  placed  in  a  tube  loop  within  the 
chamber.  The  fuel  was  forced  through  the  injector  by  opening  a 
solenoid  valve  to  a  capacitance  chamber,  maintained  at  a  pressure 
greater  than  the  test  chamber  pressure.  The  test  terminated  when 
all  the  fuel  in  the  loop  had  passed  through  the  injector.  The  in¬ 
jector  was  a  Spraying  Systems  Company  straight  hole  injector  with  a 
model  0.000009  solid  stream  tip.  The  spray  was  ignited  and  stabilized 
by  a  hot  wire  passing  close  to  the  injector  exit.  The  combustion 
process  was  observed  through  windows  in  the  test  chamber.  Motion 
picture  shadowgraphs  were  obtained  in  order  to  measure  the  boundaries 
of  the  spray.  Dark-field  motion  pictures  were  obtained  for  flame 
shape  measurements.  The  high  pressure  test  conditions  are  summarized 
in  Table  8. 

Theory.  A  preliminary  version  of  the  LHF  model  for  combusting 
sprays  was  based  on  an  integral  model  of  the  flow  and  low  pressure 
thermodynamics  (16,46).  This  earlier  model  has  been  superceded  by 
the  k-e-g  turbulence  model,  using  high  pressure  thermodynamics,  and 
will  not  be  discussed  here. 

The  basic  turbulence  model  remains  unchanged  from  the  version 
considered  in  Section  4.2.  Consideration  of  combusting  sprays  only 
requires  the  construction  of  a  new  equation  of  state.  Under  the  LHF 
assumption,  complete  thermodynamic  equilibrium  is  specified  for  each 
point  in  the  flow.  At  low  mixture  fractions,  where  the  ideal  gas 
equation  of  state  is  valid,  equilibrium  was  determined  from  the 
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Table  7 

Summary  of  Test  Conditions  for  Combusting  Jets 
at  Atmospheric  Pressure4 


Fuel 

n-Propane 

n-Pentane^ 

Fuel  State 

vapor 

liquid 

Fuel  flow  rate  (g/s) 

0.18 

0.35 

Air  flow  rate  (g/s) 

- 

0.080 

Injector  pressures  (kPa) 

Gas  side 

131.0 

120.6 

Liquid  side 

- 

'  137.9 

Injector  Thrust  (mN) 

15.6 

5.00 

Jet  Velocity  (m/s) 

88.7 

11.6 

aAmbient  conditions  296  K,  97  kPa,  injector  diameter  1.194  mm, 
injector  oriented  vertically  upward. 

^Sauter  mean  diameter  35  pm  measured  at  x/d  =  170,  at  the  centerline, 
in  cold  flow. 
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Table  8 

Summary  of  Test  Conditions  for  Combusting 
Sprays  in  High  Pressure  Air 


Fuel 

n-pentane 

Sauter  mean  diameter 

28  uma 

Maximum  drop  diameter 

79  yma 

Spray  velocity 

68. 5  m/s 

Ambient  temperature 

20  C 

Ambient  pressure 

3,  6  and  9  MPa 

aCalculated,  cf  Ref.  16 
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NASA  CEC  72  program  (83).  At  high  mixture  fractions  and  high  pressures, 
compressibility  and  solubility  effects  are  important  and  were  treated 
using  the  Redlich-Kwong  equation  of  state  (57).  For  hydrocarbons, 
some  decision  must  be  made  concerning  soot  formation  at  high  mixture 
fractions.  This  was  handled  by  assuming  that  no  soot  was  formed  at 
temperatures  below  1000  K. 

An  illustration  of  a  typical  equation  of  state  for  a  combusting 
spray  appears  in  Fig.  29.  The  system  involves  the  combustion  of  an 
n-pentane  spray,  from  a  pressure  atomizing  injector,  in  air.  This 
pressure  is  6  MPa  which  is  roughly  1.8  times  the  critical  pressure  of 
pure  n-pentane.  High  pressure  effects  are  evident  at  large  mixture 
fractions.  In  the  region  where  liquid  is  present,  f  <0.8,  the  vari¬ 
ation  of  all  species  concentrations  is  nonlinear.  This  is  due  to  the 
fact  that  significant  quantities  of  gas  are  dissolved  in  the  liquid 
phase.  Near  the  stoichiometric  mixture  ratio,  f  =  0.066,  gas  tempera¬ 
tures  reach  a  maximum,  the  density  is  a  minimum,  and  the  concentration 
of  oxygen  and  fuel  approach  zero.  For  the  complete  range  of  mixture 
fractions,  the  density  varies  almost  two  orders  of  magnitude. 

Results  and  Discussion.  Throughout  all  calculations,  the  empirical 
parameters  used  in  the  model  were  not  changed  from  the  values  given  in 
Table  6.  Only  the  equation  of  state  changes  from  case  to  case. 

The  mean  velocities  and  temperatures  measured  along  the  centerline 
of  the  combusting  gas  jet  are  illustrated  in  Fig.  30.  The  predictions 
of  the  theory  are  also  presented  on  the  plot.  The  agreement  between 
theory  and  experiment  is  excellent.  The  comparison  was  similar  for 
profiles  of  mean  velocity  and  temperature  in  the  radial  direction. 

Based  on  these  results,  it  is  clear  that  the  model  is  capable  of  handling 
combustion  processes  in  gas  flows,  with  no  additional  modification  of 
constants. 

The  comparison  between  predictions  and  measurements  for  the  com¬ 
busting  n-pentane  spray  at  atmospheric  pressure  is  illustrated  in 
Fig.  31.  Mean  velocities  and  temperatures  are  plotted  as  a  function 
of  distance  from  the  injector.  In  this  case,  there  is  evidence  of  the 
model  overestimating  the  rate  of  development  of  the  flow.  There  was 
not  sufficient  time  available  during  the  current  investigation  to 
complete  calibrated  drop-like-histories  for  this  case.  However,  the 
injector  was  the  same  as  for  the  evaporating  spray,  and  operating 
conditions  are  similar,  drop  sizes,  length  of  process,  rate  of  decay  of 
mean  velocity,  etc.  Therefore,  effects  of  slip,  like  those  pictured 
in  Fig.  27,  would  be  expected  for  the  combusting  spray  illustrated  in 
Fig.  31.  Similar  to  the  case  of  an  evaporating  spray,  the  radial 
variation  of  mean  velocity  and  temperature  was  predicted  reasonably 
well. 


The  comparison  between  the  predictions  of  the  model  and  spray 
penetration  length  measurements  for  the  high  pressure  combusting 
pentane  spray  are  summarized  in  Table  9.  The  degree  of  agreement 
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30  Theoretical  and  experimental  mean  velocity  and  temperature 
variation  along  the  centerline  of  a  combusting  gas  jet  at 
atmospheric  pressure. 
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x/d 

Fig.  31  Theoretical  and  experimental  mean  velocity  and  temperature 
variation  along  the  centerline  of  a  combusting  spray  at 
atmospheric  pressure. 
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Table  9 

Comparison  of  Measured  and  Predicted  Spray 
Penetration  Lengths  at  High  Pressuresa 


Pressure  (Mpa) 

Spray  Length/ Injector  Diameter 

Measured 

Predicted 

3 

132 

128 

6 

119 

79 

9 

100 

67 

^Test  conditions  appear  in  Table  8. 
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between  predictions  and  measurements  is  encouragingly  good.  However, 
spray  penetration  length  is  not  as  sensitive  a  measure  as  mean 
temperature  and  velocity  comparisons.  Therefore,  further  experi¬ 
mental  results  are  needed  to  provide  a  convincing  demonstration  of 
the  application  of  the  LHF  model  at  high  pressures.  Drop-life-history 
calculations  should  also  be  completed,  in  order  to  more  rigorously 
assess  the  effects  of  slip  and  loss  of  chemical  equilibrium  between 
the  phases. 

Conclusions .  The  model  was  very  effective  in  predicting  mean 
velocities  and  temperatures  in  a  combusting  gas  jet.  An  important 
feature  of  this  comparison  is  that  the  constants  used  in  the  model  are 
the  same  as  those  generally  in  use  for  turbulent  jet  computations  (76). 
These  constants  were  evaluated  again  in  the  present  investigation  by 
systematic  comparison  with  measurements  for  constant  density  and 
variable  density  jets.  The  fact  that  no  additional  parameter  was 
adjusted  to  obtain  the  comparison  shown  in  Fig.  30,  other  than  a 
standard  chemical  v,-quilibrium  calculation  to  obtain  an  equation  of 
state,  suggests  that  the  model  has  some  potential  for  a  priori  pre¬ 
dictions. 

There  was  evidence  of  slip  effects  for  the  atmospheric  pressure 
spray,  although  the  analysis  did  provide  qualitative  agreement  with 
the  measurements.  The  predictions  of  spray  penetration  lengths  for 
high  pressure  combustion  were  reasonably  good,  although  this 
parameter  is  a  relatively  course  measure  of  model  accuracy.  Further 
experimental  results  and  the  completion  of  drop -life-history  calcu¬ 
lations  are  needed  for  a  more  thorough  assessment  of  the  LHF  model 
for  combusting  sprays. 

The  advantage  of  the  approach  outlined  here  is  that  ther i  is  a 
minimum  of  "art"  required  to  complete  a  prediction.  Only  the  injector 
geometry,  flow  rate,  momentum  and  the  initial  state  of  the  reactants 
must  be  specified.  An  equation  of  state  must  be  constructed,  using 
conventional  equilibrium  thermodynamics.  All  other  parameters  of  the 
solution  are  fixed.  The  solution  yields  a  qualitatively  correct 
picture  of  the  combustion  process,  at  least  in  the  cases  considered 
to  date.  However,  the  solution  tends  to  overestimate  the  rate  of 
development  of  the  flow,  since  finite  rate  processes  between  the 
phases  are  ignored.  Thus,  the  results  of  the  LHF  prediction  represent 
a  lower  bound,  all  real  sprays  will  develop  more  slowly.  Completing 
drop- life-history  calculations,  using  the  LHF  model  results  to 
estimate  local  conditions  within  the  spray,  provides  a  systematic 
means  of  evaluating  the  potential  accuracy  of  the  LHF  model  for  a 
given  application. 
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